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ABSTRACT 

Following on from a recently completed radial- velocity survey of the old (7 Gyr) open cluster NGC 
188 in which we study in detail the solar-type hard binaries and blue stragglers of the cluster, here we 
investigate the dynamical evolution of NGC 188 through a sophisticated A-body model. Importantly 
we employ the observed binary properties of the young (150 Myr) open cluster M35, where possible, 
to guide our choices for parameters of the initial binary population. We apply pre-main-sequence 
tidal circularization and a substantial increase to the main-sequence tidal circularization rate, both 
of which are necessary to match the observed tidal circularization periods in the literature, includ- 
ing that of NGC 188. At 7 Gyr the main-sequence solar-type hard-binary population in the model 
matches that of NGC 188 in both binary frequency and distributions of orbital parameters. This 
agreement between the model and observations is in a large part due to the similarities between the 
NGC 188 and M35 solar-type binaries. Indeed, among the 7 Gyr main-sequence binaries in the model, 
only those with P > 1000 days begin to show potentially observable evidence for modifications by 
dynamical encounters, even after 7 Gyr of evolution within the star cluster. This emphasizes the 
importance of defining accurate initial conditions for star cluster models, which we propose is best 
accomplished through comparisons with observations of young open clusters like M35. Furthermore, 
this finding suggests that observations of the present-day binaries in even old open clusters can provide 
valuable information on their primordial binary populations. However, despite the model's success 
at matching the observed solar-type main-sequence population, the model underproduces blue strag- 
glers and produces an overabundance of long-period circular main-sequence - white-dwarf binaries 
as compared to the true cluster. We explore several potential solutions to the paucity of blue strag- 
glers and conclude that the model dramatically underproduces blue stragglers through mass-transfer 
processes. We suggest that common-envelope evolution may have been incorrectly imposed on the 
progenitors of the spurious long-period circular main-sequence - white-dwarf binaries, which perhaps 
instead should have gone through stable mass transfer to create blue stragglers, thereby bringing both 
the number and binary frequency of the blue straggler population in the model into agreement with 
the true blue stragglers in NGC 188. Thus, improvements in the physics of mass transfer and common 
envelope evolution employed in the model may in fact solve both discrepancies with the observations. 
This project highlights the unique accessibility of open clusters to both comprehensive observational 
surveys and full-scale A^-body simulations, both of which have only recently matured sufficiently to 
enable such a project, and underscores the importance of open clusters to the study of star cluster 
dynamics. 

(galaxy:) open clusters and associations: individual (NGC 188) - (stars:) binaries: spectroscopic - 
(stars:) blue stragglers - (methods:) numerical 



1. INTRODUCTION 

NGC 188 is one of the most well studied open 
clusters, and is a corne rstone of the WIYN Open 
Cluster Study fWOCS: iMathieul l2000|). Its old 
age (~7 Gyr; ISaraiedini et al.l I1999L iMeibom et al.l 
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I 2009D. rich binary population (|Geller et al.l 120091: 
iGeller fclvl athicu 2012) and abundance of blue straggler s 
(BSs: iMathieu fc Gellerl [200l [Geller fc Mathieul [2011 . 
X-ray sources dBelloni et al.l 119981: iGondoinl 120051 . and 
other exotic stars (see e.g.. lGeller et al E00al2009l) make 
NGC 188 particularly important for the study of star 
cluster evolution. 

Sophisticated A^-body codes now incorporate stellar 
dynamic s and stellar evolution self-consistently (e.g., 
NB0DY6; lAarsethl 120031 ). enabling models of rich open 
clusters, like NGC 188, that include up to a few 10 4 single 
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and multiple stars initia lly (e.g.. iPortegies Zwart et al"] 
1200 It IHurlev et ai1l2005| ). Until recently, the observa- 
tions needed to test many of the important parameters 
in these models, and particularly those of the binary 
population, have been unavailable, due in part to the 
decade(s) of multi-epoch observations necessary for de- 
tecting and defining the binary stars. Our comprehen- 
sive radial- ve locity survey of the sing l e and binary stars 
in N GC 188 (IGeller et al.ll2lM [200a IMathieu fc Gelled 
2009t lGeller fc Mathieull2011ll2012h now provide this es- 
sential empirical guidance, permitting us to model NGC 
188 with a high level of accuracy. 

Perhaps the most detailed iV-body study of a specific 
open cluster to date is that of lHurlev et~all (|2005t) . They 
aimed to model the open cluster M67, and particularly 
to match the observed color-magnitude diagram (CMD), 
with a focus on the BS population. Hurley et al. ini- 
tialized their binary population so as to maximize the 
production rate of BSs, and also to attempt to recreate 
the general binary properties of the M67 BSs. In doing 
so, they managed to nearly match the observed number 
of BSs in M67 at 4 Gyr in their simulation, and they 
reproduced the true cluster CMD quite closely. 

As both NGC 188 and M67 are old open clusters (at 
~7 and ~4 Gyr, respectively) and the binary and BS 
populations of bo t h clus ters are remarkably similar, in 
IGeller fc Mathieul (^Oll ) we compare this simulation to 
the binary and BS populations observed in NGC 188. 
We find that the M67 simulation contains far too many 
short-period binaries as compared to our observations of 
NGC 188, and that the binary orbital parameters and 
hard-binary (periods < 10 4 days) frequency of the BSs 
in the simulation do not match those of NGC 188 (or 
M67). These discrepancies are in a large part the result of 
the unrealistic initial binary population that was heavily 
weighted towards short-period systems, inconsistent with 
observed binary popula tions in you ng open clusters, like 
M35 ( e.g. lMathieull2008l and see also lMeibom fc Mathieul 
f2005h . 

In this paper, we model NGC 188 with a primary goal 
of matching the observed solar-type binary population. 
We employ the observed binary population of the young 
open cluster M35 (180 Myr) to guide our definition of the 
initial binary frequency and distributions of period and 
eccentricity. Thus here we create the first ./V-body model 
whose initial binary population is defined by such de- 
tailed observations of a young open cl uster. We first dis- 
cussed the BSs created in this model in lGeller fc Mathieul 
(|2011l) . Here we present the NGC 188 model in full de- 
tail. 

In Section [2] we discuss the simulation method. Sec- 
tion |31 provides a brief outline of the observations used 
in guiding the model, and Section [4] explains the initial 
conditions of the model, including a detailed discussion 
of the initial binary population. In Section [5] we analyze 
the 7 Gyr cluster CMD, mass and structure in compar- 
ison to observations of NGC 188. We then compare the 
binary population and BSs in the model at 7 Gyr to 
the observed binaries and BSs in NGC 188 (Sections |5] 
and [7]) . In Section [8] we analyze the dynamically formed 
triples within our model as compared to observations of 
triples in open clusters and the Galactic field. Finally in 
Section [9] we investigate the dynamical evolution of the 
main-sequence (MS) binaries and evaluate BS produc- 



tion within our A-body model, and we then provide our 
conclusions in Section [TOl 

2. SIMULATION METHOD 

We use the NBDDY6 code (|Aarsethll20"03[ ) to model the 
dynamical evolution within the NGC 188 simulations. 
Stellar and binar y evolution a r e incl uded in NBODY6 
using the work of lHurlev et al.l (|2000[ ) and IHurlev et al.l 
$2&m. Much of our method is identical to that of 
IHurlev et al.l ((2005), with one distinction being that they 
used the NB0DY4 code wi th GRAPE-6 computing hard- 
ware (Ma kino et al"1l2003ft while we use the NBDDY6 code 
on the supercomputer at the Swinburne Centre for Astro- 
physics and Supercomputing (without a GRAPE board). 
The main difference between the two code versions is 
that NBDDY4 interfaces with the GRAPE hardware for 
the force calculations while in NBDDY6 this is done on a 
standard CPU or more recently on a gra phics process- 
ing unit fGPU; iNitadori fc Aarsethl 12012ft . To increase 
efficiency, NB0DY6 also has the option to compute force 
contributions from di stant particles less frequently than 
from near neighbors (| Aarsethl Il999l . We modified the 
NB0DY6 code to define the initial binary population and 
output format. We choose to output stellar evolution 
parameters, binary orbital properties, positions and ve- 
locities for all stars in snapshot intervals of ~30 Myr. 

Individual stars are evolved according to the Single- 
Star Evolution algorithm fSSE; IHurlev et al.ll2000D . This 
code rapidly models all phases of evolution from the zero- 
age main sequence (ZAMS) through remnant phases, 
for stars covering a mass range of 0.1 - 100 M Q and 
a wide range in metallicity. Mass loss through stellar 
winds is included. Binary stars are evolved according to 
the B inary-Star Evolution algorithm (BSE: lHurlev et all 
2002). This code contains prescriptions for tidal cir- 
cularization and synchronization, angular momentum 
loss mechanisms (e.g., magnetic braking and gravita- 
tional radiation) , mass transfer from Roche lobe overflow 
(RLOF), accretion from stellar winds, common-envelope 
(CE) events, mergers, and direct collisions. Thus the 
code contains numerous pathways to create anomalous 
stars, such as BSs, even without the aid of dynamical 
encounters within a star cluster. The additional consid- 
erations that may occur in a cluster environment, such as 
perturbations t o binary orbits , exchanges, etc., are mod- 
el ed in NB0DY6 (lAarseth 2003, and references therein) . 

IHurlev et all |2002j, |2005|) explain the modeling of BSs 
in detail. As a goal of this paper is to closely compare 
the observed and simulated BS populations, we provide a 
short summary here. BSs can form in the A-body model 
through three primary mechanisms, namely mass trans- 
fer during RLOF, physical stellar collisions, and mergers 
of two stars within a binary. BS formation through wind 
mass transfer is also possible in principle, but in practice 
this formation channel is very infrequent. Therefore for 
this paper, we will use the term mass transfer to refer 
specifically to mass transfer from RLOF. 

Mass tra nsfer during RLOF is modeled in NB0DY6 ac- 
cording to IHurlev et al.l (|2002l ). In short, the stability 
of mass transfer is det ermined u s ing th e radius-mass ex- 
ponents £ defined by IWebbinkl (|1985f ). which are used 
to distinguish between mass transfer that proceeds on a 
dynamical, thermal or nuclear time scale. BSs created 
through mass transfer from RLOF derive from thermal, 
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and sometimes nuclear mass transfer, and are formed 
when a MS star accretes mass from either a MS (Case 
A), red giant (Case B) or asymptotic giant (Case C) 
donor. Case A mass transfer leads to a merger, while 
Cases B and C generally result in a BS in a detached 
binary with a white dwarf (WD) companion. Dynamical 
mass transfer either leads to a CE episode, when mass 
is transferred from a giant, or a merger, when mass is 
transferred from a low-mass (< 0.7 M Q ) MS star, nei- 
ther of which will produce a BS at 7 Gyr. In prac- 
tice the criteria to determine whether a binary under- 
going RLOF will enter dynamical mass transfer or ther- 
mal/nuclear mass transfer (e.g., possibly creating a BS) 
is determined by the critical mass ratio, q c , such that if 
Mdonor/M accret or > 1c the binary undergo es dynamical 
mass trans fer. Both the q r formu l ae fro m iHurlev et al.l 
(|2002f) and iHiellming fe Webbinkl (|1987[ ) are options in 
the A-body code, and we have chosen to implement the 
latter. We return to the discussion of CE evolution and 
the q c value in Section 19721 

If a MS star gains mass through a mass-transfer 
episode, merger or collision, it is presumed to remain on 
the MS but attain a higher luminosity and effective tem- 
perature. The lifetime of the star is determined based on 
the fraction of unburned hydrogen in the core of the star, 
such that the fraction of MS lifetime remaining is directly 
proportional to the fraction of remaining unburned hy- 
drogen in the core. 

During thermal and nuclear mass transfer, the rejuve- 
nation procedure is determined based on the structure 
of the accretor star's core. If mass is transferred onto a 
star with a radiative core (0.35 M < M < 1.25 M ) the 
fraction of unburned hydrogen is barely changed. There- 
fore the fraction of its MS lifetime that has elapsed is 
unaffected by the mass transfer, but the effective age of 
the star is less than that of other normal stars of simi- 
lar mass. If mass is transferred onto a star with a con- 
vective core (M > 1.25 M Q ) or that is fully convective 
(M < 0.35 Mq), the core grows and mixes in unburned 
hydrogen. It is assumed that the amount of burnt hydro- 
gen is unchanged, and that the core mass grows directly 
proportional to the mass of the star, which is increased 
during mass transfer. The increase in the fraction of 
unburned hydrogen in the core increases the remaining 
fraction of MS lifetime. If mass transfer causes a star 
to move across the radiative/convective-core boundary 
the treatment simply switches from one core type to the 
other. After the mass-transfer phase, the total MS life- 
time of this rejuvenated star (from t = in the simulation 
to the time when the star evolves off the MS) for either 
stars with convective or radiative cores is equal to the 
total lifetime of a normal star of that mass. The rejuve- 
nation process described here ensures that the fraction of 
MS lifetime remaining for this star is greater than that 
of a normal MS star of this mass (regardless of whether 
the core is radiative or convective). 

If two MS stars undergo a merger or collision, the prod- 
uct is taken to be a new MS star where the stellar ma- 
terial is fully mixed. A BS formed in this manner then 
evolves along the MS similarly to normal MS stars of 
the same mass until 10% of the total hydrogen has been 
burnt, at which point the BS evolves off of the MS and 
follows the standard evoluti onary sequen c e for a star of 
this mass. As discussed in IHurlev et al.l ()2005l ) the as- 



sumption of complete mixing here is likely incorrect and 
may affect the lifetime and appearance of these stars. We 
follow this method here in part to aid comparisons with 
previous models that use similar prescriptions. Though 
correcting for this over simplification is highly desired, it 
is beyond the scope of this project. We also note that, in 
the A-body code, collision and merger products are as- 
sumed to achieve thermal equilibrium rapidly, and these 
BSs follow equilibrium models upon creation. However, 
detailed studies of collision products show that such BSs 
may not be in thermal equilibrium, and may be more lu- 
minous and have short er lifetimes than the e quilibrium 
models used here (e.g. Glcbbec k et al.l [2008). We will 
return to this point briefly in Section [3J 

The processes involved in forming a BS can cause the 
star to lie above the MS turnoff on a CMD. Observation- 
ally we identify BSs as being generally both brighter and 
bluer than the MS turnoff. Therefore here we will also 
use this definition, and not include any collision, merger, 
or mass-transfer products that still remain below the MS 
turnoff. This empirical identification procedure neces- 
sarily introduces some ambiguity between binaries found 
near the turnoff, as a normal equal-mass MS-MS binary 
could be up to 0.75 magnitudes brighter than the turnoff 
(without containing a BS). Within the simulation we can 
also identify BSs as being more massive than a normal 
MS star at the turnoff, while still maintaining the struc- 
ture and evolutionary state of a MS star. Therefore this 
ambiguity does not effect BS identification in the A-body 
model. 

In order to reduce the stochastic effects present in N- 
body simulations, and especially those of BS production, 
we ran twenty simulations, and for much of the analysis 
we average the results together. For all twenty simula- 
tions, the initial stellar and binary parameters (e.g., po- 
sitions, velocities, binary periods and eccentricities, etc.) 
are all drawn from the same respective distributions, but 
each simulation randomizes the parameters to produce a 
unique initial stellar population. In this paper we will use 
the term "simulation" to refer to any individual A-body 
simulation (e.g., any one of these twenty realizations), 
and we will use the term NGC 188 "model" to refer to 
the combination of these twenty simulations. We dis- 
cuss our method for combining these simulations in later 
sections. 

3. OBSERVED BINARY STAR POPULATIONS 

Here we give a brief outline of the observed stellar 
samples in M35 used to define our initial binary pop- 
ulation and in NGC 188 for comparison with the 7 
Gyr binary and BS populations in the A-body model. 
T he M35 radia l -veloc ity sample is discussed in detail 
in iGeller et al.l (|2010|) . In short, our M35 sample in- 
cludes MS solar-type stars within the magnitude range of 
13.0 < V < 16.5 (1.6 - 0.8 M©) and within 30 arcminutes 
from the cluster center (7 pc in projection at a distance 
of 805 pc or ~4 core radii). The cluster turnoff in M35 is 
at V ~9.5 (~4 M(7)), and the r efore our sample includes 
only MS stars. IGeller et all (|2010l ) identify 360 single 
and 55 binary cluster members within this sample; 39 of 
these binaries have orbital solutions. We use this binary 
population to define the initial binary frequency and dis- 
tributions of binary orbital parameters for the NGC 188 
simulations (Section I4.2|) . 
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The NGC 188 rad i al-vel ocity sample is discussed in 
detail in IGeller et a l. (2008). Our stellar sample spans a 
magnitude range of 10.5 < V < 16.5 (1.1 - 0.9 M Q ), and 
extends to a projected radius of 17 pc (^13 core radii, or 
30 arcminutes). The sample includes the brightest stars 
in the cluster down to about two magnitudes below the 
MS turnoff, and all NGC 188 BSs. We identify 358 sin- 
gle and 129 binary cluster members within this sample; 
85 of these binaries have orbital solutions. The distri- 
butions of orbital param eters are discussed in detail in 
IGeller fe Mathieul (|2012t ) and we also show a sample of 
these results in Section [Bl 

For both clusters we can detect binaries out to peri- 
ods o f order 10 4 days. The hard-soft boundary (Hcggic 
1975) in both of these clusters is at periods of order 
~ 10 6 days, and therefore all detected binaries in these 
clusters are hard binaries. Here we will use the term 
"hard binary" to refer specifically to binaries within the 
period range that we can detect with our observations. 
The observed binary frequencies and distributions of or- 
bital parameters for M35 and NGC 188 presented here 
ha ve been corrected for inco mpleteness using the method 
of IGeller fc Mathieul (|2012t ). and importantly, our com- 
pleteness in binary detection and orbital solutions is 
roughly equivalent for both clusters. 

4. INITIAL CONDITIONS OF THE SIMULATIONS 
4.1. Initial Cluster Structure 

NGC 188 is obs erved to have a solar m etallicity and 
an age of 7 Gvr (ISaraiedini et all Il999l and see also 
iMeibom et"al"1l2009l and iFornal et alJl2007| ). Empirical 
estimates for the mass of NGC 188 range from about 
1500 M^ to 3800 Ma, (IBonatto et alj |2005t IGeller et alJ 
120081 : IChumak et al II2010D . In order to ret ain this mass 
by 7 G yr, we scale from the results of the IHurlev et al.1 
( 2005) simulation to account for mass loss through evapo- 
ration and stellar evolution processes. We find a suitable 
initial mass of 23600 Mq (3 9000 stars). W e choose the 
stellar masses according to a iKroupal (|2001l ) initial-mass 
function (IMF). 

We distribute the stars in space according to a Plum- 
mer density profile. In theory the Plummer model ex- 
tends to infinite radius; therefore in practice we apply 
a cut-off at ten times the half-mass radius to avoid rare 
cases of stars placed at large distances in the initial dis- 
tribution. This initial Plummer density profile quickly 
evolves to resemble a King profil e within our model (as 
in other jV-body models, e.g. IHurlev fc Sharal [20031 : 
IHurlev et ail 12005). The initial velocities are isotropic 
and are generated based on the positions and the as- 
sumption of dynamical equilibrium. 

The initial size scale of a simulation is defined by the 
initial half-mass radius, and we aim to match the ob- 
served core and tida l radii of the true cluster at 7 Gyr. 
IBonatto et all (|2005| ) find NGC 188 to have a tidal ra- 
dius of 21 ±4 p c and a core radius of 1.3±0.1 pc, while 
IChumak et "all (|2010h find a tidal radius of MA±{ 6 pc 

and a core radi us of 2.llg | pc (both using standard 
KingfllllQjll model fits). We note that IChumak et all 
(2010) included only proper- motion membe rs of th e clus- 
ter in their analysis while IBonatto et al.1 (|2005f l used 
2MASS photometry without any kinematic membership 
data, which may account for the discrepancies between 



these two results. For our purposes, we assume that the 
true core and tidal radii of NGC 188 fall somewhere be- 
tween these values. 

We ran a number of test simulations with different half- 
mass radii to attempt to estimate an appropriate initial 
half-mass radius that will evolve to match the observed 
values described above. We find that using an initial half- 
mass radius of 4.6 pc results in a reasonable agreement 
with the observed core and tidal radii for NGC 188 at 
the age of 7 Gyr, see Section [SJ and we use this initial 
half- mass radius f or ou r NGC 188 model. (Also note 
that IHurlev et al.1 (|2005f ) used an initial half- mass radius 
of 3.9 pc, and found a resulting core radius of 0.64 pc at 
4 Gyr, much smalle r than NGC 188.) 

ICarraro fc Chiosil (|l99l find NGC 188 to be in a 
nearly circular orbit (e = 0.07) at a dis tance between 
9.4 to 10.8 kpc from the Galactic center. IChumak et al.1 
(|2010l) find the cluster to orbit between 6.97 and 9.03 
kpc with a mean Galactocentric distance that is nearly 
the same as that of the Sun. They also find the or- 
bit to have a low eccentricity of e — 0.125. The differ- 
ences between these two results are due to the Galactic 
models used and som e updated observations used in the 
IChumak et al.l (|2010f) analysis. 

For simplicity we place the cluster in a standard 
Galactic tidal field for the Solar neighborhood, having 
a circular orbit with a speed of 220 km s _1 at a dis- 
tance of 8.5 kpc from the Galactic cente r. Adding the 
non-ze ro orbital eccentricity suggested bv lChumak et ahl 
(2010) would introduce periods of increased stellar es- 
cape rate due to the reduced tidal radius during peri- 
galacticon, as well as the opposite effect at apogalacticon. 
iMadrid et al.1 (|2012| ) show that, after a Hubble time, a 
cluster in an elliptical orbit with an eccentricity of 0.5 and 
apogalacticon of 8 kpc will have ~10% less mass than the 
same cluster in a circular orbit at 8 kpc from the Galac- 
tic center. Thus w e assume that includi ng the small ec- 
centricity found bv lChumak et al.l (|2010T ) would result in 
only a minimal change to the mass of the model at 7 Gyr. 
We also note that NGC 188 is currently out of the plane 
of our Galaxy , at b = +22.°38, an d its orbit crosses the 
galactic disk (|Chumak et al.ll2010D . We do not account 
for any effects related to this disk crossing (e.g., tidal disk 
shocks) in our model. Disk shocks temporarily compress 
a cluster on a timescale shorter than the internal dy- 
namical timescale, which accelerates all stars, some to 
velocities above the escape v elocity (mostly those in the 
outer parts of the cluster). IVesperini fc Heggiel ()1997f ) 
show that for a cluster orbiting at 8.5 kpc, disk shocks 
produce only a mini mal effect on the clus ter mass. This 
result is confirmed bv lKiipper et al.l (|2010l ) who also show 
that the lifetime of a cluster on an 8.5 kpc orbit is much 
more strongly effected by a highly eccentric orbit than a 
highly inclined orbit. 

4.2. Parameters of the Initial Binary Population 

The initial parameters defining a binary population in 
numerical simulations have long been derived from the- 
oretical analyses and predictions, as we have lacked suf- 
ficiently detailed observations of young open cluster bi- 
nary populatio ns to define these pa rameters. The WOCS 
survey of M35 (jGeller et al.l feOlO') now allows for a ma- 
jor step forward; we will use these binaries to empirically 
guide our choices for the initial binary population in these 



Direct TV-body Modeling of NGC 188 



5 





Fig. 1. — Initial distributions of solar-type binary orbital parameters for the NGC 188 model. The distributions of period (top-left), 
eccentricity (top-right), secondary mass (bottom-left) and mass ratio (bottom-right) are shown in histogram form, and the eccentricity 
- log period (e — log(P)) diagram is shown in the center. Solar-type binaries within our observational limits (primary mass from 0.8 - 
1.6 Mq and P < 3000 days) from the model are shown in the gray (green in the online version) filled histograms and points. In the 
period and e — log(P) distribution plots we also show the full solar-type (0.8 - 1.6 M©) sample in the thin-lined histogram and small 
points, respectively. The o bservable sample in the eccentricity distribution plot excludes binaries with P < P c i rc , where P c i rc = 10.2 days 
(Mcibom & Mathicu 2005). Each histogram bin shows the mean value for all simulations included in the model, and the error bars show 
one standard deviation above and below the mean. The e — log(P) figure shows one representative simulation. The theoretical functions 
from which we draw the initial binary periods and eccentricities in the model are shown in the black curves in the upper two plots (i.e., 
the log-normal period distribution from Duqucnnoy & Mayor 1991 an d the Gaussian fit to the M35 eccentricity distribution). We describe 
our method for choosing the initial mass-ratio distribution in Section 14.21 Finally the observed M35 period and eccentricity distributions 
are plotted in thick-black histograms. 



simulations. 

We don't yet know if all clusters, including NGC 188, 
form with similar binary populations. However, as we 
will show in Section [6j an initial binary population based 
on observations of M35 will indeed match our observa- 
tions of the NGC 188 binaries at 7 Gyr. We discuss the 
implications of this result in Section [§1 

The M35 binaries in our sample show a rising distri- 
bution in period that agrees with the log-normal distri- 



butio n observed for solar-type binaries in the Galacti c 
field (|Duauennov fc Mavonfl99lt iRaghavan et al.ll2010D . 
As we only detect binaries with periods < 10 4 days, we 
utilize these Galactic field binaries to extend the period 
distribution to longer periods. Specifically we input a 
primordial binary population with a log-normal period 
distribution c entered on log(P [days]') = 4 .8 with a = 2.3 
(as defined bv lDuauennov fc Mavorl ll991). We don't an- 
ticipate the initial soft binary population to survive long 
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against dynamical encounters, and so our results will be 
only minimally sensitive to this extrapolation. 

We observe M35 to have a MS b inary frequency o f 
24% ± 3% out to periods of 10 4 days (1 Geller et al.ll2010D . 
We choose to set the initial binary frequency in the model 
to 27% out to the same period cutoff. We allow the bi- 
naries to initially populate the entire period distribution, 
which corresponds to an initial total binary frequency of 
^70% (16000 binaries). However, the initially very long- 
period soft binaries (P > 10 days) are quickly broken 
up through dynamical interactions. Thus the total bi- 
nary frequency drops rapidly to 53% within the first ^30 
Myr (see Figure [7]) . 

The M35 binary eccentricity distribution (for binaries 
with Pdrc < P < 3000 days, where P circ = 10.2±^ 
days from iMeibom fc Mathieul [2005) is consistent with 
a Gaussian distribution centered on e = 0.38 with 
a = 0.23. Thus we draw the initial binary eccentrici- 
ties from this Gaussian function. We note that this ec- 
centricity distribution is again similar to that observed 
for solar-type binaries in the field, xs We have not 
yet performed an analysis of the observed secondary- 
mass or mass-ratio distributions for the M35 binaries. 
These distributions are generally the least certain ob- 
servationally, at least for spectroscopic samples, as the 
majority of observed spectroscopic binaries are single 
lined. Observations of the Galactic field binary popu- 
lations show the mass-ratio distribut ion to be either ris- 
ing to wards lower ma ss ratios (e.g. IDuauenno v fc Mayor! 
19911) or uniform (e.g. lMazeh et a l. 2003; Rag havan et al.l 
20101) . There is also ev idence for a peak in the mass-ratio 



distri b ution at unity (iRaghavan et al.l 120101 : iTokovininl 
12000ft . iReggiani fc Meverl (|201lD find a mass-ratio distri- 
bution in their binary sample that follows a power law 
dN/dq oc q&, with /3 = -0.50 ± 0.29, for binaries with 
primary masses between 0.25 M Q and 6.5 M Q . 

We choose to define the binary mass es by fir s t tak - 
ing two masses, Mi and M 2 , from the iKroupal ()2001f ) 
IMF, both more massive than 0.1 Mq. This is approxi- 
mately the hydrogen-burning limit, and no models below 
0.1 Mq were evolved for the SSE code. We then com- 
bine these masses (M to t) and randomly choose a new 
mass ratio from a uniform distribution to define the pri- 
mary and secondary masses for the binary, such that 
Mi + AI2 = M to t and the secondary star has a mass 
M 2 in the range Mi < M 2 < 0.1 M Q . The result of this 
procedure is shown in the bottom panels of Figure [1] 

This method does not produce the observed peak at 
mass ratios of unity for solar-type binaries; we will inves- 
tigate here whether this peak can be produced through 
dynamical processes throughout the evolution of the clus- 
ter. The solar- type binary mass-ra tio distribution in our 
model agrees with the results of IDuauennov fc Mayor! 
(|199lD . who studied solar-t ype binaries in th e field , 
and also with the results o f Rcggiani fc Meverl (|2011[ ) 
and lMetchev fc Hillenbr and (20Q9|) who studied a wider 
range of binaries in different environments, all of whom 
find that the distribution of mass ratios increases to- 
wards lower mass rati os. (Although we note that 
IReggiani fc Meverl 1 2 1 ll do not find a variation in the 
mass-ratio distribution between different mass samples, 
which does result from our pairing method.) We there- 
fore proceed with this method for defining the component 



masses of the binaries in the NGC 188 model, but we 
note that future models may benefit from a more thor- 
ough investigation of the effects of different methods for 
choosing binary component masses on the evolution of 
the binary population. 

The initial binary distributions are shown in Figure [U 
In this figure, as well as throughout the paper we distin- 
guish between the "observable" and "total" samples such 
that the observable sample includes only stars that are 
within the appropriate mass and period range for our 
radial- velocity observations of the true cluster(s), and 
total refers to the entire simulated cluster population. 
Here the observable distributions are shown in filled his- 
tograms and larger points. 

4.2.1. Tidal Energy Dissipation Rate 

We digress briefly here to discuss the tidal treatment 
in the iV-body code and the additional parameters that 
we impose on the model to match our observations. The 
tidal circularization period is one of our best observa- 
tional to ols to study the effec t s of ti des on a binary pop- 
ulation. IMeibom fc Mathieul (|2005l ) define the circular- 
ization period as the period at which the best-fit circular - 
ization function (equation 1 in lMeibom fc Mathieull2005l) 
reaches an eccentricity of 0.01. They use this method to 
study the circularization period of a number of binary 
populations of different ages, and we have reproduced 
the majority of their Figure 9 here in Figure [5] with the 
labeled points. They find a trend of increasing circular- 
ization period with increasing age of the binary popula- 
tion, for_ages_>l Gyr. 

InKJellcr et al. (2009) we note that the tidal treatment 



in the TV-body model produces circularization periods 
that are well below those that are observed in real open 
clusters. Indeed using the standard tidal energy dissipa- 
tion rates in the iV-body code results in a circularization 
period for solar-type binaries at 7 Gyr of ~4 days, as 
compared to the obs e rved v alue of 14.5+2;2 days found by 
IMeibom fc Mathieul ()2005D for similar binaries in NGC 
188. Therefore here we investigate methods for improv- 
ing the agreement between the circularization periods 
derived from the iV-body model and those from the ob- 
served binary populations throug h a series of popula tion 
synthesis simulations using BSE dHurlev et alJl2002l ). 

The solid line in Figure [2] shows the circularization pe- 
riod as a function of time for a sample of isolated binaries 
evolved using the standard tidal prescription from the N- 
body and BSE codes. We derive all circularization peri- 
ods from the po pulation synthesis simulati ons using the 
same method as IMeibom fc Mathieul (|2005l ). and then fit 
exponential functions to these discreet circularization pe- 
riods to produce the smooth curves shown here. The cir- 
cularization periods from the standard BSE and A-body 
codes fall below those of nearly all observed clusters. 
These low circularization periods suggest that the tidal 
energy dissipation rate may be under estimated within 
the m odel. This has also been noted bvlBelczvnsk i et all 
(2008), who suggest that the convective tidal damping 
rate should be increased by a factor of 50. We find that 
a factor of 100 increase is necessary in our model to be 
consistent with the lower limit on the tidal circularization 
period of NGC 188 at 7 Gyr (dashed line). 

However for ages less than ~4 Gyr the circularization 
period is still quite low as compared to observations of 
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Fig. 2. — Circularization period a s a function of age. The la - 
beled points reproduce the data from Mcibom & Mathicu (2005), 
and we also add in the lower limit for the circular i zation period 
of NGC 6819 (2.4 Gyr) found by IMilliman et al.l pOlI) , The 
lines show fits of exponential functions to the results from mul- 
tiple p opulation synthesis studies using the binary-star evolution 
(BSE; IHurlev et al.l [20021) algorithm. Specifically, the solid line 
shows the default BSE settings. The dashed line shows the result 
of increasing the convective damping term by a factor of 100. The 
dotted line shows the result of also add i ng th e pre-main-sequence 
circularization algorithm from [Kroupa (1993). The default BSE 
settings result in circularization periods that are consistently less 
than observations of real binary population at any age. With the 
addition of the pre-main-sequence circularization algorithm as well 
as the increased convective damping strength, the model matches 
the observed circularization periods. 

M35, the Pleiades and the pre -MS binary population. In- 
deed, IZahn fc Bouchet] (|1989f ) predict that pre-MS tidal 
circularization will result in circularization periods be- 
tween 7.2 and 8.5 days, in agreemen t with the empirical 
results of IMeibom fc Mathieul (|2005l ). 

In order to achieve this early tidal circularization, we 
ad d the pre-MS b inary evolution prescription suggested 
by IKroupal (|1995l ) . 1 The results of this combination of 
both the increased conv ective tidal damping coefficient 
and the IKroupal (|1995h pre-MS binary evolution pre- 
scription is shown in the dotted line. (For comparison, 
using the pre-MS binary evolution prescription without 
increasing the convective tidal strength results in essen- 
tially the same P C irc ~ 8 days at all times.) 

We note that this procedure is simply the ad hoc ad- 
dition of two mechanisms, and here we use them to- 
gether only to (a) achieve the desired initial form to 
the e — log(P) distribution and (b) maintain a cir- 
cularizatio n period that is con s istent with the obser- 
vations of IMeibom fc Mathieul (|2005| ) out to the age 
of NGC 188. Future work is necessary to investigate 
whether such a "hybrid" style model can be developed 
into a more com prehensive tidal the ory (as has also 
been suggested by j Mathicu et al.|[i"992l and discussed by 
IMeibom fc Mathiei]ll20r)5r i~ 

In summary, to initialize the binaries in our NGC 188 
model, we first choose orbital parameters independently 

1 In order to better fit the "envelope" in observed e — log(P) 
dist ribution s of op en clusters, we modify the A and \ values in 



the IKroupal (|1995l 1 algorit hm from 28 and 0.75, to 20 and 15, re- 
spectively. As described in Kroupa (1995), A measures the length 
scale over which significant evolution of the orbital elements dur- 
ing the proto-stellar phase occurs, and \ measures the "interaction 
strength" between the t wo protostars in the binary system. See 
IGeller fc Mathieu (2012) for a comp arison of the e — log(P) dis- 
tribution in the IHurlev et al.l (120051 ) model resulting from the A 
and x values from IKroupa ( 1995J with the observed e — log(P) 
distribution of NGC 188. 



from the observationally denned distributions discussed 
above. In rare cases we modify these initial orbital pa- 
rameters to avoid unphysical systems, (e.g. those that 
would collide or merge in a single orbit). We then mod- 
ify theljhianesaccording to the pre-MS binary evolution 
of IKroupal (|1995h . The resulting distributions of orbital 
parameters are shown in Figure Q] As the binaries evolve 
within the model they are subject to the increased tidal 
damping factor described above. 

5. COLOR-MAGNITUDE DIAGRAM, MASS AND 
STRUCTURE OF THE NGC 188 MODEL AND 
OBSERVATIONS AT 7 GYR 

We evolve the stellar population, as defined above, dy- 
namically with the NBDDY6 code to an age of 7 Gyr to 
create a simulation of NGC 188. As discussed above, we 
ran twenty simulations each with a different realization 
of the same initial parameter distributions, and for much 
of the analysis that we discuss here and in the following 
sections, we combine the results of these twenty simula- 
tions to form the NGC 188 model. 

First we discuss the 7 Gyr CMD, Figure |3l In the 
left panel of this figure, we show one representative sim- 
ulation; the other 19 simulations have qualitatively the 
same features, although the specific numbers of stars in 
various populations (e.g. BSs and triples) vary. Here we 
briefly compare this simulation to the observed CMD for 
NGC 188, shown in the right panel of same figure. 

The magnitudes for each star in a given simulation ar e 
calculated in the same manner as in IHurlev et al1([2005T l. 
using the bolometric co rrections o f Kuruc zl (I1992D for MS 
and giant stars and of iBergeron et al.l (|1995l ) for WDs. 
We use the distance modulus (m — M)y = 11 .44 and red - 
dening E(B-V) = 0.09 found bv lSarajedini et al.l (|1999l ) 
to place the simulation at the appropriate distance and 
color. In both CMDs, we show all stars with the small 
black points and highlight binaries within our observa- 
tional detection limits with larger points, using light gray 
for MS binaries, dark gray for red giant branch (RGB) 
binaries, and black for BS binaries (green, red and blue in 
the online version). We use this same gray-scale (color) 
coding in our analyses of these respective binary popu- 
lations throughout the paper. 

At the distance of NGC 188, the components of a circu- 
lar binary at the hard-soft boundary (P ~ 10 6 5 days for 
solar-type stars at 7 Gyr) would be separated by < 0.01 
arcseconds on the sky. Therefore, no binaries (or higher 
order systems) are resolvable from ground-based obser- 
vations, and here we treat all binaries (and higher order 
systems) in the model as unresolved objects. 

The model qualitatively reproduces the MS single and 
binary sequences quite closely as compared to the ob- 
served CMD. We see a width to the MS in both CMDs 
as a result of the rich binary population of the cluster. 
(We will discuss the binaries in detail in Section [6l) The 
cluster turnoff is at V ~ 14.8 in both CMDs, which cor- 
responds to a mass of ~1.1 Mq. The MS luminosity 
function from the model is also consistent with that of 
the proper- motion members from lPlatais et ail (|2003| ) in 
both number and form (for V < 19, below which the 
iPlatais et al.ll200"3l catalog becomes incomplete). 

However, there are two obvious discrepancies between 
the observed and simulated clusters. First the model 
does not predict the large long-known scatter on the 
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Fig. 3. — Color-magnitude diagrams (CMDs) of a representative simulation from the NGC 188 model (left) and observations of the 
true cluster (right). For both panels, the detectable binaries (V < 16.5 and P < 10 4 days) are plotted in the larger filled circles with the 
main-sequence binaries in light gray, giant binaries in dark gray and blue straggler binaries in black (green, red and blue in the online 
version). Solar- type binaries with periods longer than 10 4 days in the simulation are plotted in open circles. Triples are plotte d in crosses. 
All other stars are plotted in black points. In the observed CMD, we show all cluster member s from the Geller et al. (2008) survey and 
extend the sample to fainter magnitudes using the proper-motion members (Ppm > 50%) from Platais ct al. (2003). The proper motions 
extend to a faint limit of V = 21, but begin to become incomplete at V ~ 19. For the simulation, binaries and triples are plotted as single 
points, showing t heir c ombined magnitudes and colors. The distance modulus of (m — M)v = 11-44 and reddening of E(B — V) = 0.09 of 
ISaraiedini et af] 1 119991 ) are employed to place the simulation at the same distance and color as NGC 188. 

Next we investigate the cluster structure at 7 Gyr. For 
this and later analyses we analyze each simulation in the 
NGC 188 model as it might appear to an observer on 
Earth. NGC 188 has Galactic coordin ates of I = 122.°85 
and b = +22.°38 (jPlatais et al.ll2003l ). Each simulation 
begins with the X-direction facing away from the Galac- 
tic center and the cluster orbits in the X-Y plane. We 
rotate each simulation according to the Galactic coordi- 
nates such that the transformed Y' — Z' plane becomes 
the plane of the sky, and the X' direction is the line-of 
sight. We then use this observed line-of-sight to NGC 188 
in order to investigate the cluster structure at 7 Gyr, and 
for subsequent analyses of the model where appropriate. 

In Figure [5] we show the stellar surface density as a 
function of radius from the cluster center at the age of 
7 Gyr. Binaries and higher-order systems are treated 
as individual objects. Each point in this figure shows 
the mean stellar density at the given radius, averaged 
over all simulations. The error bars show the standard 
deviations of these values . Th i s prof ile is fit with the 
empirical three-parameter [King (1962;) model, as shown 
by the solid black line (xl e d = 1-9)- The fit results in a 
core radius of 1.52±0.08 pc and a tidal radius of 25±3 pc. 



[19771 iTwarod [19781 IGeller et all [20081) . The source of 
this observed scatter is not due to the precision of the 
photometry. The BV photometr y shown here is precise 
to 0.01 mag (jPlatais et al.ll2003f ). We have not added 
photometric errors to the simulated CMD. The origin of 
this observed scatter on the N GC 188 RGB, al so seen in 
other clusters (e.g., NGC 6819, lHole et al.ll2009D . remains 
unknown. 

Second, the model clearly has fewer BSs at 7 Gyr than 
the true cluster. We will discuss the BSs in detail in 
Section [7] 

Next, in Figure 0] we show the mean total mass of 
the model as a function of time. The mean initial mass 
of all twenty simulations is 23610 M e , with a standard 
deviation of 40 Mq. At 7 Gyr the simulated clusters 
have lost on average ^88% of this mass through evap- 
oration processes and stellar evolution, leaving a mean 
final mass of 2850 Mq, with a standard deviation of 
120 Mq. This value is in good agreement with em- 
pirical mass estimates for NGC 188 that range from 
1500 to 3800 Mq (Igonatto et alJl2005b [Geller et al.l 
l2008t IChumak et al.ll20ll )~ 
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Fig. 4. — Logarithm of the mean total cluster mass as a 
function of time in the NGC 188 model. The black line is cen- 
tered on the mean total mass at a given time from our twenty 
simulations, and the width of the line shows one standard de- 
viation above and below the mean. The light-gray horizontal 
band shows the range in mass derived e mpirically f or NGC 188 
of 1500 M(p to 3800 M Q (IBonatto et al.ll2005l ; [Geller et alj|2008t 
IChumak et al. 20j3). The model begins with a mean total mass 
of 23610 Mq (with a standard deviation of 40 Mq between the 
simulations), and at 7 Gyr the mean total mass is reduced to 
2850 Mq (with a standard deviation of 120 Mq) through stellar 
evolution and evaporation processes. The mean mass of the model 
at the age of NGC 188 agrees with the observed mass estimates. 

With the dotted and dashed black lines in Figure [SJ 
we show two recent results analyzing the observed stel- 
lar density pr ofile in NGC 188. Th e dotted line shows 
the results of iBonatto et all (|2005| ) who used 2MASS 
data to derive a core radius of 1.3 ± 1 pc and a tidal 
ra dius of 21 ± 4 p c. Th e dashed line shows the results 
of IChumak eTail $2Q M who use the NGC 188 proper- 
motion members from iPlatais et al.1 (|2003| ) to derive a 
core radius of 2.1+gg pc and a tidal radius of 34.4j^Q pc. 
The values for the NGC 188 iV-body model fall comfort- 
ably between the results of these two studies. 

In Figure [6] we compare the cumulative radial distri- 
butions of the observable single, binary, RGB and BS 
popu lations in th e mod el with those of the true cluster 
from lGeller et all (|2008t ). We include the BS radial pro- 
file here for completeness, but we will wait to discuss the 
BS population until Section [3 For Figure [6] and all sub- 
sequent cumulative distribution plots presented in this 
paper, we use the union of all twenty simulations to con- 
struct the cumulative distribution functions for the NGC 
188 model. The solar-type binaries in the model are cen- 
trally concentrated with respect to the solar-type single 
stars at the >99% confidence level, in agreement with 
the true cluster. This result also holds for any individual 
simulation. The solar-type binaries are on average more 
massive than the single stars, and this result reflects the 
fact that the model is mass segregated at 7 Gyr, as has 
also been observed in the true cluster (ISaraiedini et al.l 
[19991 iKafka fe Honevcuttl[2Q0l [G^Iier et al.ll2008lh ~ 

The RGB stars in the model appear to be more cen- 
trally concentrated than the single stars when combining 
all twenty simulations (as shown in Figure [5]). This is in 
contrast to the observed radial distribution of the RGB 
stars in NGC 188, which follow a nearly identical radial 
distribution to the single stars in the true cluster. How- 
ever for any given simulation, a K-S test does not return 
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Fig. 5. — Stellar surface density as a function of radius from the 
cluster center at 7 Gyr in the model, projected into the true line- 
of-sight towards NGC 188. Each point shows the mean value for 
all simulations within the given radial bin, and the error bars show 
one standard deviation above and below the mean. We include 
only stars (and binaries) from the simulation with V < 21, to 
match the limiting magnitude of the Platais et al. (2003) proper- 
motion study of the true cluster. The solid black line shows a 
standard King model fit to these data, which yields a core radius 
r c = 1.52 ± 0.08 pc and a tidal radius of rt = 25 ± 3 pc. For 
comparison we also show two King model fits from the literature 
to the observed dens ity profile of NGC 188. The King model fit of 
IBonatto et a l. (2005) is shown in the dotted line, who find NGC 
18 8 to have a r c = 1.3 p c and rt = 21 pc, and the King model fit 
of IChum ak ct al. (2010) is shown with the dashed line, who find 
r c = 2.1 pc and rt = 34.4 pc (both normalized to the same central 
density as our fit to the NGC 188 model). The vertical line shows 
the mean half-mass radius of the model at 7 Gyr. 

a significant distinction between the radial distributions 
of the RGB and single stars (given the relatively small 
number of RGB stars per simulation). 

6. THE 7 GYR MAIN-SEQUENCE AND GIANT BINARY 
POPULATIONS 

6.1. Binary Frequency 

At 7 Gyr in the twenty simulations the mean observ- 
able MS hard-binary (P < 10 4 days) frequency is 33.5%, 
with a standard deviation of 2.8%. This result is in 
agreement with the observed MS hard-binary frequency 
of NGC 188 of 29% ± 3% (within th e same period and 
mass range; [Geller fc Mathieul [201 2) . As is clear from 
Figure [7] the model predicts that the binary frequency 
will change only slightly over 7 Gyr of evolution. More- 
over, a cluster with a binary frequency consistent with 
that of M35 at 180 Myr will evolve to have a binary fre- 
quency consistent with NGC 188 at 7 Gyr. 

The mean hard-binary frequency of the RGB stars in 
the model is 31% with a standard deviation of 6%, in- 
distinguishable from that of the model MS stars. This 
bin ary frequency also agrees with that of the true clus- 
ter; [Geiir&Sathie^ |20ll) find the NGC 188 giants to 
have a hard-binary frequency of 34% ± 9%. 

Again we see clear evidence of mass segregation (in 
both the simulated and true cluster) when we examine 
the hard-binary frequency of the observable MS and RGB 
populations combined as a function of radius from the 
cluster center (Figure [8|). We exclude the BSs from the 
analysis, because the radial dependence of their binary 
frequency may be tied to their formation mechanisms. 
Here the model is show n in the filled circles and is com- 
pared to the result from lGeller fc Mathieul (|2012f) for the 
true cluster MS and RGB stars, shown in open circles. 
The agreement is remarkable. 
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Fig. 6. — Radial cumulative distri butions of cluster p opulations from the NGC 188 model (left; in observed line-of-sight projection) and 
observations (right; reproduced from Geller et al. 2008|). The cluster is divided into the observable singles (solid line), observable binaries 
(dotted line), giants (dashed line) and blue stragglers (dot-dashed line). The solar- type binaries are centrally concentrated with respect to 
the solar-type single stars in the model and observations, due to two-body relaxation processes. The simulated blue stragglers do not show 
the same bimodal spatial distribution that is observed in NGC 188. 
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Fig. 7. — Binary frequency as a function of time in the NGC 188 
model. The total binary frequency is shown with the black line, 
the total main-sequence binary frequency with the light-gray line, 
and the observable main-sequence binary frequency with the dark- 
gray line (M> 0.9 Mq and P < 10 4 days). All lines are centered 
on the mean values of all twenty simulations, and the widths are 
one standard deviation above and below each mean. Also shown 
(as points) are the observed main-sequence hard-binary frequen- 
cies of M35 (24% ± 3% at 180 Myr) and NGC 188 (29% ± 3% 
at 7 Gyr), both for P < 10 4 day binaries. The observable main- 
sequence binary frequency in the model is consistent with those of 
both clusters at their respective ages. We also note that the total 
binary frequency drops below the main-sequence binary frequency 
as the cluster evolves due to the population of binaries with white- 
dwarf primaries, which has a low binary frequency. (In detached 
binaries, mass-loss on the giant branch increases the orbital separa- 
tion and thereby increases the probability for dynamical encounters 
resulting in exchanges or disruption.) 

6.2. Distributions of Periods, Eccentricities, Secondary 
Masses and Mass Ratios 

In Figures 1^1 - IT51 we show the eccentricity - log period 
(e — log(P)) distribution and the distributions of eccen- 
tricities, periods, secondary masses and mass ratios for 
the observable solar-type binaries in the model as com- 
pared to observations of similar binaries in NGC 188. In 
all figures, the MS, RGB and BS binaries from the model 
are plotted in light gray, dark gray and black (green, red 
and blue in the online version) . For the selection of the 
RGB and BS binary populations, only one star must be 
a RGB or BS star, respectively. For completeness, we in- 
clude the BSs in Figuresl!?l-[T3l but we wait to discuss the 
BS population and the selection of the BS sample until 
Section [7] Here we focus on the MS and RGB binaries. 



Fig. 8. — Binary frequency as a function of radius at 7 Gyr for 
the NGC 188 model and observations. For both the model and 
observations, the sample includes main-sequence and giant stars 
combined. Filled circles show the model binary frequencies and 
open circles show the observed binary frequencies. Errors bars 
on the observations show the Poisson counting uncertainties. The 
points from the model show the mean values from all simulations, 
and the error bars show one standard deviation above and below 
the means, respectively. The horizontal bars indicate the range 
in radius from the cluster center for which the binary frequency is 
calculated. (The bins used for the observed and simulated data are 
identical; the plotted points are shifted along the x-axis for clar- 
ity.) The binary frequency increases towards the core in the model 
due to two-body relaxation processes. Remarkably, the binary fre- 
quency in the model matches that of the true cluster at each bin 
in radius. 

In Figure |9] we show the e — log(P) diagram for the 
observable MS (top), RGB (middle) and BS (bottom) 
stellar populations in the model. The MS and RGB 
e — log(P) distributions a re in close agreement wi th the 
observations (Figure 3 in iGeller fc Mathieu! I20T2I) . The 
mean 7 Gyr circularization period over all of our sim- 
ulations is 17.3 ± 2.1 da ys, consistent with t h e valu e of 
U.5±H days found by 'iMeibom fc Mathieu! (f2005h for 
NGC 188. 

However, we notice a significant population of cir- 
cular MS binaries with P > P c irc in the model, 
which is also seen in the eccentricity distribution shown 
in Figure 1101 Such binaries ar e not observed in 
NGC 188 (IGeller fc Matlfieul [jOll), the Ga l actic field 
(|Duauennov fc Mayor! 119911: iRaghavan et ail I2010D . or 
any of the open cluster b inary populations studied by 
IMeibom fc Mathieu! (|2005| ). We will return to these bi- 
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Fig. 9. — Eccentricity plotted against the logarithm of the pe- 
riod for the main- sequence (top), giant (middle) and blue strag- 
gler (bottom) binaries in the model. For the main-sequence and 
giant binaries, we show all binaries present in all simulations at 7 
Gyr. The integrated sample of blue stragglers shown in the bottom 
panel contains all blue-straggler binaries present in the model at 
each snapshot interval be tween 6 and 7.5 Gyr. The circularization 
period P c irc = 14.5 days JMeibom fc M athicu 2005j) is plotted for 
reference with the dashed vertical line in the main-sequence panel. 
In all three panels, we circle inner binaries of triple systems, and 
plot the outer orbital parameters for triples with crosses. 



naries in Section [9.21 Here we simply note that on aver- 
age in a given simulation < 1 of these long-period circular 
binaries had an initial eccentricity < 0.01. The vast ma- 
jority started with non-zero eccentricities, but were circu- 
larized during a CE episode which left a WD secondary 
star. If we do not include these post-CE systems, the 
MS eccentricity distribution of the model is statistically 
indistinguishable from the observed distribution. 

At 7 Gyr the solar-type MS binary period distribu- 
tion in the model is also consistent with that observed 
in NGC 188. The simulated MS period distribution in 
Figure [Til rises towards our observational detection limit 
as does the observed distribution. Including only bina- 
ries that are within our observational completeness limit 
of P < 3000 days, the observed and simulated period 
distributions are statistically indistinguishable. 

The MS secondary-mass and mass-ratio distributions 
in the model (Figures [12] and [T3l) have forms some- 
what different from those of the true cluster. Both dis- 
tributions in the true cluster rise towards lo wer-mass 
companions (although iGeller fc Mathieull2012l note that 
they cannot formally distinguish the observed mass- 
function distribution, the observable quantity, from 
a distribution derived by choosing binary component 
masses from a uniform mass-ratio distribution). In 
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Fig. 10. — Eccentricity distributions for the main-sequence, gi- 
ant and blue strag gler binaries with P c j rc < P < 3000 days (where 
P circ = 14.5 days; IMeibom fc Mathicu 200^). The period limit is 
set empirically to exclude binaries whose orbits have been circu- 
larized by tides, and to extend to our completeness limit in orbital 
solutions for the NGC 188 sample. In the top three panels the 
model is shown in filled histograms, showing all solar- type main- 
sequence and giant binaries present at 7 Gyr and the integrated 
blue straggler sample. Each histogram bin represents the mean 
number of binaries in all twenty simulations within the given ec- 
centricity range, and the error bars show the standard errors of 
these means. For comparison we plot the observed NGC 188 dis- 
tributions in the thick-lined histograms, respectively. The observed 
NGC 188 blue straggler frequency distribution is divided by 5 for 
clarity, and the true observed frequency is shown on the right y- 
axis. In the bottom panel we plot the cumulative distributions for 
the main- sequence, giant and blue straggler samples in the model 
in light gray, dark gray and black (green, red and blue in the online 
version) . For comparison we also plot the solar- type main-sequence 
eccentricity distribution of the initial population (thin solid line) 
and that resulting from evolving the same initial binaries in isola- 
tion for 7 Gyr (dashed line). 

both the MS secondary-mass and mass-ratio distribu- 
tions from the model, we see distributions that rise until 
M 2 ^0.5 M Q and q ~0.5, followed by a turnover towards 
lower masses and mass ratios. The peak in both distri- 
butions is in part a result of WD companions, the vast 
majority of which are members of the unphysical long- 
period circular MS- WD binaries discussed briefly above. 
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Fig. 11. — Period distributions for the main- sequence, giant 
and blue straggler binaries. The plots are of the same format as 
Figure [10] except here we do not limit the model by period. (The 
observations extend to a period of 3000 days.) Additionally we 
show the full log-normal period distribution from which we chose 
the initial binary periods in the top panel, normalized to the 7 Gyr 
solar-type main-sequence binary frequency. 



If we remove these MS- WD binaries (dotted lines in the 
MS panels of Figures [T^] and [T3J) , the distributions in 
secondary-mass and mass-ratio are both consistent with 
a uniform distribution, respectively. 

We see a similar peak at about 0.5 M in the 
secondary-mass and mass-ratio distributions for the 
RGB binaries. This peak is also primarily due to WDs. 
Furthermore ~60% of these RGB - WD binaries evolved 
from post-CE MS - WD binaries with long-periods and 
circular orbits, similar to those discussed above. Nearly 
all of the remaining 40% are also evolved post-CE MS 
-WD binaries, but with P < P C i rc . Again, we will dis- 
cuss these long-period circular post-CE MS- WD binaries 
in detail in Section |9~21 

The observed NGC 188 MS mass-ratio distribution 
shows a peak at mass-ratios near unity (also seen in 
the higher-mass MS stars in the observed secondary-mass 
distribution), commonly referred to as "twins" and also 




M 2 ( M , 

Fig. 12. — Secondary-mass distributions for the main- sequence, 
giant and blue straggler binaries. The plots are of the same for- 
mat as Figure [TOl except here we include all solar- type observable 
binaries with P < 3000 days. Additionally in the main-sequence 
panel we show the distribution excluding the unphysical MS- WD 
binaries with the dotted line. Note the plots begin at 0.1 Mq, as 
this is t he minimum ini tial stellar mass evolved within the TV-body 
code. (Geller & Mathieu (2012) do not derive a secondary-mass 
distribution for the giants.) 



observed in the Galactic field (e.g. lRaghavan et aLll2010| ). 
The model does not create such equal mass ratio systems, 
which suggests that these binaries are likely not formed 
dynamically and instead are a result of the binary for- 
mation process. 

Recall that we have not yet derived the secondary- 
mass and mass-ratio distributions for the solar-type stars 
in M35. Therefore these initial distributions are uncer- 
tain in our model. Investigating the cumulative distri- 
butions in Figures [12] and [13] shows that the observ- 
able secondary-mass and mass-ratio distributions in the 
model have only changed slightly from their respective 
initial conditions. These results suggest that NGC 188 
may have formed with a secondary-mass distribution 
that was more strongly weighted towards lower masses 
(possibly following the IMF directly), and also a mass- 
ratio distribution with a peak at mass ratios near unity. 
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Fig. 13. — Mass-ratio distributions for the main-sequence, giant 
and blue straggler binaries. The plots are of the same format as 
Figure H"2l 

Finally, we note that the period, eccentricity and 
secondary-mass distributions of the RGB binaries in the 
model are all statistically indistinguishable from the re- 
spective MS distributions. The model is consistent with 
the observed RGB period and eccentricity distributions, 
although we note that the sm all observed sample size p ro- 
hibits a precise comparison. (IGeller fc MathieullMa did 
not derive the observed RGB secondary-mass or mass- 
ratio distributions for NGC 188.) 

6.3. The Role of Stellar Dynamics in Shaping the 
Main-Sequence Binary Population 

About 70% of the MS-MS binaries present in the model 
at 7 Gyr maintained their primordial periods and eccen- 
tricities (both to within 10%), despite 7 Gyr of dynamical 
evolution within the star cluster. Furthermore, 84% of 
the 7 Gyr solar-type MS-MS binaries within our observ- 
able mass and period limits (P < 3000 days, M\ > 0.9 
M Q ) maintained their primordial periods and eccentric- 
ities. In short, the distributions of orbital parameters 
for the MS binaries, and especially for short-period bi- 
naries, were largely maintained throughout the cluster 
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Fig. 14. — Mean eccentricity as a function of orbital period for 
solar-type main-sequence - main-sequence binaries in the model 
(with initial binaries in light gray and 7 Gyr binaries in dark gray) 
and observations (black). We show binaries with P > P c i rc , and 
plot the mean eccentricity over the bins in period shown in the 
horizontal lines. Vertical error bars for the observed data show the 
standard errors of the mean. For the model, we plot the mean val- 
ues across all twenty simulations, and the vertical error bars show 
one standard deviation above and below these means. At periods 
> 1000 days, the 7 Gyr mean eccentricities become significantly 
larger than the initial eccentricities, due to stellar encounters. A 
hint at a similar increase in eccentricity with increasing period is 
seen in the NGC 188 observations. 

evolution. 

In order to further investigate how dynamics has mod- 
ified the binary population over the 7 Gyr of cluster evo- 
lution, we compare the initial and 7 Gyr solar-type MS 
period, eccentricity, secondary-mass and mass-ratio dis- 
tributions with those derived by evolving the binaries 
from the model in isolation for 7 Gyr. In the cumulative 
distribution panels of Figures ll0l - [T3l the initial distribu- 
tions are shown in the solid lines. The 7 Gyr solar-type 
MS distributions evolved within the cluster are shown 
in the light-gray (green in the online version) points, and 
the 7 Gyr MS distributions evolved in isolation are shown 
in the dashed lines. 

Given the union of twenty simulations, we are able to 
distinguish the initial and 7 Gyr MS binary period, ec- 
centricity secondary-mass and mass-ratio distributions, 
respectively, all at very high confidence. Likewise we 
can distinguish at very high confidence these 7 Gyr bi- 
nary distributions from the respective distributions de- 
rived from isolated evolution. However considering any 
individual simulation, only the 7 Gyr orbital period dis- 
tribution has changed significantly from the initial distri- 
bution, with a significance level of >99%. Similarly the 
7 Gyr MS period distribution evolved within the clus- 
ter can be distinguished from that evolved in isolation at 
the >99% confidence level for any given realization of the 
cluster. All similar comparisons for the other binary dis- 
tributions for any single simulation result in significance 
levels well below 90%. 

Only the longest period binaries experience significant 
orbital modifications from dynamical encounters. In the 
case of the period distribution, we began our model with 
a fully populated log-normal distribution (solid line Fig- 
ure Hip . Dynamical encounters have the most dramatic 
effect on the longest period, soft binaries, which are 
quickly disrupted by dynamical encounters, imposing a 
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hard-soft boundary at P ~ 10 6 5 days. The result is a pe- 
riod distribution at 7 Gyr that can be distinguished from 
the initial distribution at high confidence. Conversely, 
the binaries evolved in isolation maintain the full log- 
normal period distribution and show no distinction from 
the initial distribution. 

In the case of the eccentricity distribution, we began 
with a Gaussian distribution at all periods. Again, we see 
that dynamical encounters affect the longest period bina- 
ries most severely (Figure [TJ]). The eccentricity distribu- 
tions for binaries with periods > 1000 days are shifted to 
significantly higher eccentricities than their initial dis- 
tributions, and the difference increases with increasing 
period. Wider binaries have larger cross sections for en- 
counters, and such encounters (most often fly-by encoun- 
ters) tend to increase eccentricity. We note, however, 
that even binaries with periods very close to the hard- 
soft boundary do not attain a thermal eccentricity dis- 
tribution in these simulations. 

In Figure [14] we also compare to observations of the 
NGC 188 MS binaries (black points). We see a hint 
of a similar trend towards higher eccentricities amongst 
longer period binaries in the observations. However 
given the relatively small sample size, this possible trend 
in the observations is n ot significant. As discussed in 
IGeller fc Mathieul (l20ll . many of these higher eccen- 
tricity binaries in NGC 188 reside near the cluster core, 
where dynamical encounters are most frequent. 

Stellar dynamical and relaxation processes have also 
resulted in a modest increase in binary frequency (Fig- 
ure [7]) . The mean 7 Gyr solar- type MS hard-binary fre- 
quency over all simulations is 33.5% (with a standard 
deviation of 2.8%), an increase of about 20% over the ini- 
tial MS solar-type hard-binary frequency. This increase 
in binary frequency is due to the preferential evaporation 
of the lower-mass single stars as compared to the higher- 
total-mass binary stars in our observable mass range. 

As is clear from this analysis, the agreement between 
the 7 Gyr solar-type MS hard-binary population in the 
model with our observations is largely due to our choice 
of initial conditions, as most hard-binaries within our 
observable limit in orbital period have not been dramat- 
ically affected by strong dynamical encounters over the 7 
Gyr of evolution. We return to this point in Section 19.11 

7. COMPARISON BETWEEN THE OBSERVED AND 
SIMULATED BLUE STRAGGLERS 

For the BS population in the model, we use an in- 
tegrated sample that includes all BSs present at each 
^30 Myr snapshot interval within the twenty simulations 
from 6 to 7.5 Gyr, covering the r ange in the observe d age 
estimates for NGC 188 (e.g. see lFornal et al.ll2007ft . and 
we will do similarly for the majority the paper. The 
use of multiple simulations reduces the stochastic fluc- 
tuations in BS production. The inclusion of all BSs at 
each snapshot interval weights the BS population by the 
lifetime that a given BS spends in each cluster radius, 
orbital configuration, etc., and essentially provides the 
likelihood that we will observe a BS in a given location 
o r with a given orbi t al co nfiguration, etc. 

IGeller fc Mathieul (|20ll use the BS population from 
this model to examine the predicted binary orbital pa- 
rameters of BSs formed through mass transfer and stel- 
lar collisions and to investigate the origins of the long- 



period BS binaries in NGC 188. They conclude, based on 
the long-periods, low eccentricities and particularly the 
secondary-mass distribution which has a narrow peak at 
~0.5M Q , that the NGC 188 BSs in long-period bina- 
ries, which comprise the majority of the BS population 
in the cluster, likely have origins through mass-transfer 
processes. Here we expand upon their results. 

First we examine the radial distribution of the BSs 
(Figure The observed NGC 188 BSs show a bimodal 
radial distribution with a centrally concentrated group 
of 14 BSs near th e core and an addi ti onal population of 
7 BSs in the halo (IGeller et al.ll2008lL IGeller fe Mathieul 
((Ml) find no significant difference between the binary 
properties of these two BS population in NGC 188, sug- 
gesting that this bimodal distribution is not a result of 
two distinct BS formation channels. Similar bimodal 
BS radial distri butions have been observed in globu- 
lar clusters (e.g. iFerraro et al.lll997t ) and have been at- 
tribute d to dynamica l fricti on and mass-segregation pro- 
cesses (|Mapelli et al.l [2004]) . These same processes act 
within our NGC 188 model; however the BS population 
in the model does not show the bimodal spatial distribu- 
tion that we observe in the true cluster. Moreover only 
17% of the simulated BSs are found outside of 5 pc (~3.5 
core radii), as compared to ^33% in NGC 188. Later we 
suggest that this difference is a result of incorrect model- 
ing of BS formation channels (and particularly the mass- 
transfer channel) which does not adequately population 
the halo with BSs. 

In Figure ED we show the 21 BSs observed in NGC 188 
at 7 Gyr, 16 of which ar e in binaries with P < 10 4 days 
(|Mathieu fc Gellerll2009l) . At 7 Gyr, the maximum num- 
ber of BSs produced in any of the twenty simulations is 
11, and the mean number of BSs at 7 Gyr is 6.2 ± 0.8. In 
Figure [T5l we show the number of BSs (top) and the BS 
hard-binary frequency (P < 10 4 bottom) as functions 
of time for the NGC 188 model. We find on average 
6.09 ± 0.05 BSs between 6 and 7.5 Gyr in our integrated 
BS sample from the model, less than one third of the 
observed number of BSs in NGC 188. 

Next we compare the observed and simulated BS bi- 
nary properties, starting with binary frequency. The 
mean BS hard-binary frequency found in the model is 
also below the observed value of 76% ± 19%. Between 
6 and 7.5 Gyr, the mean detectable BS binary frequency 
is 14.5% ± 0.5%, a factor of 5 lower than the BS hard- 
binary frequency observed in NGC 188. As is clear from 
Figure Qj)] both the total number of BSs and the BS 
hard-binary frequency from the model are well below the 
observed values for essentially the entire lifetime of the 
modeled cluster. Furthermore, the entire BS frequency, 
independent of period, between 6 and 7.5 Gyr is only 
27.0% ± 0.4%, about one third of the observed value 
(which only included BS binaries with P < 10 4 days). 

We now return to Figures [10] - 1131 where we compare 
the distributions of the orbital parameters for the inte- 
grated sample of BSs with those of the observed BSs in 
NGC 188. Importantly, the forms of the distributions of 
orbital parameters for the detectable BS binaries in the 
model are consistent with those that we observe for the 
true BSs in NGC 188. A large fraction of the simulated 
BSs are found in low-eccentricity orbits with long peri- 
ods, and their secondary-mass distribution is peaked at 
~0.5 M Q . Respective K-S test show that the BS binary 
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Fig. 15. — Number of blue stragglers (top) and the detectable 
blue-straggler binary frequency (P < 10 4 days; bottom) as func- 
tions of time. In both panels, the black line shows the mean values 
at each snapshot interval found by averaging the results from all 
twenty simulations, and the gray region marks one standard devi- 
ation above and below the mean values. In NGC 188, there are 21 
observed blue stragglers, 16 of which are in hard binaries, resulting 
in a blue straggler hard- binary frequency of 76% ± 19%. In the 
model both the number of blue stragglers and the blue straggler 
binary frequency are significantly lower than the observed values. 

period, eccentricity secondary-mass and mass-ratio dis- 
tributions in the model can be distinguished from the as- 
sociated solar-type MS binary distributions in the model 
at the >99% confidence le vels. 

This echos the r esults o f IMathieu fe Gellerl (|2009f ) and 
iGeller fe Mathieul (|201lL l2012f ) who compared the ob- 
served BS and solar-type MS binaries in NGC 188. They 
find that the NGC 188 BSs in binaries have significantly 
longer periods than the MS binaries (at the >99% confi- 
dence level), and that amongst such long-period binaries, 
the BSs have significantly lower eccentricities (at the 98% 
confidence level). Furthermore, the companion- mass dis- 
tribution for the long-period BS binaries in NGC 188 is 
narrowly peaked near ~0.5 Mq, while the companion- 
mass distribution of solar- type MS binaries shows n o ev- 
idence for a similar peak. IGeller feTM athicu ( 201 lh sug- 
gest that these BS binary properties point to an origin 
in mass-transfer processes, where the companions would 
be ~0.5 M Q white dwarfs. 

We note here that ~60% of the BSs in binaries with 
P < 3000 days in the model derive from mass-transfer 
processes (and we expand upon this below). Thus the 
agreement between the observed and simulated distri- 
butions of BS binary parameters is likely a result of 
both populations having the same dominant BS forma- 
tion channel. 

In the model, we also find BSs in binaries at longer 
periods than our 3000 day completeness limit for ob- 
servations (see Figure fTTj) . Currently five BSs in NGC 
188 show no detectable radial-velocity variability, and 
are therefore considered "single". However these NGC 
188 BSs may prove to have companions in longer-period 
orbits (currently beyond our detection limit), analogous 
to the longer-period BS binaries in the model. 

Next, we analyze the BSs from the model according to 



formation mechanism. BSs in the NGC 188 model form 
through three primary mechanisms, namely direct stel- 
lar collisions, mergers and mass-transfer processes. Here 
we are concerned with the BS outcome. Therefore in the 
following we will include within the merger category any 
binaries that are in the process of mass transfer and will 
eventually merge, which are predominantly Case A mass- 
transfer binaries between two MS stars. Stellar collisions 
occur during dynamical encounters involving binaries (or 
higher-order systems), while mergers and mass transfer 
can occur within an isolated binary. In the NGC 188 
model the rates of mergers and mass transfer are also 
modified by dynamical encounters (as we discuss in Sec- 
tion EM- 

Recently a new BS formation path was proposed 
(|Ivanovall2008t iPerets fe Fabrvckvl |2009| ) . where the in- 
ner binary in a tr iple sy st em is driven to a merger by 
a combination of Kozai (1962) cycles and tidal fric- 
tion known as the KCTF mechanism. This mechanism 
is available in the jV-body model through the work of 
iMardling fe Aarsethl (|2001[) who added prescriptions to 
model Kozai cycles and tidal friction in triple stars. 

In general, we find that binaries and triples contribute 
significantly to BS formation in the NGC 188 model. 
The vast majority of the BSs formed in the model, re- 
gardless of the formation mechanism, had origins involv- 
ing primordial binaries. Furthermore, nearly half of the 
BSs present in the NGC 188 model between 6 and 7.5 
Gyr were found in dynamically formed hierarchical triple 
systems (in many cases with a primordial binary as the 
inner system) during the snapshot interval prior to be- 
coming a BS. The wider orbit of the tertiary increases 
the cross section for stellar encounters, and in nearly all 
of these systems, the BSs formed through collisions re- 
sulting from dynamical inter actions (as was pred icted by 
IMathieu fc Gellerl 120091 and iLeigh fe Sills! 1201 11 1. Thus 
binaries and triples are key to BS formation in the NGC 
188 model. 

However, we do not find the KCTF mechanism to form 
a significant population of BSs in the NGC 188 model, 
and therefore we will not include this mechanism in the 
following discussion of this section. Importantly, though, 
we did not include triples in the initial population of the 
model. We return to BS formation through KCTF in 
Section O 

In Figure [16] we show the number of BSs as a func- 
tion of time in the NGC 188 model, separated by for- 
mation mechanism, with BSs formed by collisions in the 
solid line, mergers in the dashed line and mass-transfer 
processes in the dotted line. BS formation early in the 
cluster lifetime is dominated by mergers of initially very 
short-period binaries. At ~3 Gyr the cluster core be- 
gins to contract, and we see a corresponding rise in the 
collision rate. At the age of NGC 188 collisions are the 
dominant BS formation mechanism, followed closely by 
mergers. The mass-transfer rate remains roughly con- 
stant throughout the lifetime of the cluster, and is in gen- 
eral significantly lower than either the collision or merger 
formation rates. 

The relative contributi ons of the dif f erent BS forma- 
tion mechanisms in the iHurlev et all (|2005l ) M67 N- 
body model are very si milar to what we find here (see 
IGeller fc Mathieul I2012D . Of the 20 BSs present at 
4 Gyr in the M67 simulation, eight formed through 
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Fig. 16. — Number of blue stragglers as a function of time in the 
NGC 188 model for blue stragglers formed through collisions (solid 
line) mergers (dashed line) and mass transfer (dotted line), binned 
into 250 Myr intervals. Each bin is normalized to show the mean 
number of blue stragglers per simulation at a given time formed 
through the respective mechanisms. Mergers dominate early in the 
cluster lifetime. At ~3 Gyr collisions begin forming BSs at a rate 
equivalent to mergers. This corresponds with the time when the 
core begins to contract in the model. The number of blue stragglers 
formed through mass transfer is always significantly lower than 
those created through the other mechanisms. At the age of NGC 
188 collisions and mergers dominate blue straggler formation in the 
TV-body model. 

collisions, ten formed through mergers, and only two 
fo rmed through mass-tran sfer processes. As discussed 
in lGeller fc Mathieul (|2012ft the M67 simulation contains 
a hard-binary frequency {P < 10 4 days) at 4 Gyr that is 
about twice that observed in NGC 188 (or included in our 
NGC 188 model), and an overabundance of short-period 
binaries. This large excess of binaries is responsible for 
the hi gher formation rate of BSs through all mechanisms 
in the iHurlev et al.l (|2005| ) simulations. 

Next we compare the binary frequencies for BSs 
(within the integrated sample) formed by each forma- 
tion mechanism in the NGC 188 model. Only 4% of the 
merger BSs in the NGC 188 model have binary com- 
panions with orbital periods below our observational de- 
tection limit. The majority of these are in contact sys- 
tems that will eventually merge to form single BSs, while 
a small fraction have exchanged into detached binaries 
with wider separations. Only 12% of the BSs formed by 
collisions in the model would be detected as binaries. In- 
terestingly, even if we consider BS binaries at any period, 
only 33% of BSs formed by collisions are found in bina- 
ries. Many of these single collisional BSs formed with 
binary companions initially in very wide orbits, near the 
hard-soft boundary, that were quickly ionized. Thus the 
low observable-binary frequency for the BSs in the NGC 
188 model reflects the dominant formation channels. 

Essentially all of the BSs formed by mass transfer in 
the NGC 188 model are in binaries that would be de- 
tectable observationally. Only the mass-transfer mech- 
anism can reproduce the large observed BS binary fre- 
quency of NGC 188. Furthermore, about 50% of the BSs 
that are detected as binaries in the NGC 188 model de- 
rive from mass-transfer processes. Thus, although the 
absolute number of BSs formed by mass transfer in the 
NGC 188 model is low, mass transfer is the most efficient 
mechanism at producing BSs in binaries. As discussed 
above, the agreement between the distributions of binary 
parameters for the BSs in the model with the BSs ob- 
served in NGC 188 is likely due to the majority of BSs in 



both samples having origins in mass-transfer processes. 

Finally, we note that the NGC 188 model struggles to 
produce BSs with similar binary characteristics as the 
two double-lined (SB2) BSs in NGC 188. Both of these 
SB2 BSs in NGC 188 have short periods (P < 10 days) 
and companions with masses Mi > 0.9 M , one of which 
is a MS star near the turnoff and one of which is itself 
a BS. There are no BS binaries with P < 10 days and 
M 2 > 0.9 M Q present in the model between 6 and 7.5 
Gyr; there are two such systems present at earlier times. 
Interestingly, out of the twenty simulations in the NGC 
188 model, we find three BS-BS binaries present between 
6 and 7.5 Gyr with short enough periods to be detected 
as binaries observationally. All of these systems formed 
through exchange encounters. However these BS-BS bi- 
naries all have periods much greater than 10 days (specif- 
ically, 340, 2750, and 5500 days, respectively), and two 
of these three systems were only bound for one snap- 
shot interval. Another 4 BS-BS binaries are present in 
the integrated sample, but with periods > 10 days, and 
therefore not detectable in our observations. To date, no 
such long-period BS-BS binaries have been detected in 
open clusters or the field. 

In summary, the deficiency in number of BSs, the low 
frequency of detectable binaries among those that are 
formed, and the lack of a bimodal BS radial distribu- 
tion are strikin g failures of the model c ompared to the 
observations of iGeller fe Mathieul (|2011l) . We return to 
this point in Section 19.31 Importantly, however, the dis- 
tributions of orbital periods, eccentricities, companion 
masses and mass ratios for the detectable BS binaries in 
the model (Figures H0]fT3|) are all consistent with the ob- 
served properties of the NGC 188 BSs. The majority of 
these detectable BS binaries in the model were formed 
through mass-transfer processes. Thus, the agreement 
between the simulated and observed distributions of BS 
bi nary orbital parameters further supports the results 
of IGeller fc Mathieul (|2011l ). who find that the majority 
of the true BSs in NGC 188 likely have origins in mass 
transfer. 

8. DYNAMICALLY FORMED TRIPLES 

Recently triple-star systems have been employed 
to explain a number of "anomalous" star systems 
observed in both star clusters and the Galac- 
tic field, including certai n bright X-ray b i naries 
(iMakarov fe Eggletonl 120091) . BSs (|Ivanoval 120081: 
iPerets fc Fabrvckvl I2009D , and some short-peri o d con - 
tact binaries (|Eggleton fe Kisseleva-Eggletonl |2006|) . 
As we did not include any triples initially, all triples 
discussed here were formed dynamically, often during 
binary-binary encounters. Observationally, the triple 
populations in star clusters (including NGC 188) are 
poorly known. Therefore here we mainly compare the 
triple frequency and distributions of orbital parameters 
in the NGC 188 model to those observed in the Galactic 
field, and investigate whether dynamical encounters 
can form triples with the same characteristics as those 
observed in the field. 

8.1. Frequency of Triples 

We first investigate the number of dynamically formed 
triples and the triple frequency (Nt r ip/N bj ec t s , where 
we treat each single star or system of stars as one ob- 
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Fig. 17. — Number of triples (top) and the triple frequency 
{Ntrip/N i,j ecta , bottom) as functions of time. We include only 
triples with solar- type primary stars. The black lines show the 
mean values at each snapshot interval for the twenty simulations, 
and the gray filled area shows one standard deviation above and 
below the means. The dotted lines show the mean values for solar- 
type triples that have outer periods < 10 4 days (our spectroscopic 
detectability limit). All triples formed dynamically, as we did not 
include triples initially in the model. 

ject, respectively) as functions of time. In Figure [T7l 
the solid black lines show the mean number (top) and 
frequency (bottom) of triples with solar-type primaries 
(0.9 < Mi(M Q )< 1.1) in the model. The dotted lines 
show the values for solar-type triples that we could po- 
tentially detect in our spectra, namely those with outer 
periods < 10 4 days. The number of solar-type triples 
increases early on due to initially very long-period bina- 
ries near (and beyond) the hard-soft boundary undergo- 
ing binary-binary encounters that result in dynamically 
formed triples, as well as an increasing encounter rate in 
the core due to core contraction. By ~5 Gyr the number 
of solar- type triples reaches a maximum of ~10 triples, 
and then declines as triples begin to evaporate from the 
cluster or are destroyed at a rate higher than they are 
created. The average lifetime of a given triple in the 
model is 820 ± 60 Myr. 

On the other hand, the solar-type triple frequency con- 
tinues to increase throughout the entire lifetime of the 
modeled cluster. These more massive systems are prefer- 
entially retained by the cluster as the less massive single 
stars are ejected. At 7 Gyr the solar- type triple frequency 
is 2.10% ± 0.13%. However many of these triples have 
o uter periods b eyond our observational detection limit. 

iGeller et al.l ([2009D report two candidate triples (and 
one quadruple system) in their sample of binaries with 
orbital solutions in NGC 188, placing a lower limit on 
the MS triple frequency of 0.5% (2/375), which certainly 
suffers from significant incompleteness. The NGC 188 
model produces a detectable solar-type triple frequency 
of only 0.15% ± 0.04% at 7 Gyr, significantly lower than 
e ven this observed lo wer li mit . 

iMermilliod et al.l (|1992l ) find a triple frequency of 
2.3% ± 1.6 % (2/88) for F5-K0 st a rs in t he Pleiades (-150 
Myr), and IMermilliod fc Mavorl (| 19991) find F5-K0 stars 
in Praesepe (~750 Myr) to have a triple frequency of 



3.8% ± 2.2% (3/80). Both of these spectroscopic surveys 
have a roughly similar completeness limit as our NGC 
188 survey. The NGC 188 model produces essentially 
zero solar-type detectable triples at the ages of these 
young clusters (as is clear from Figure [T7|). If we take 
the full sample of solar-type triples (no longer limited to 
have P < 10 4 days), we find a solar-type triple frequency 
of 0.11% ± 0.01% at 150 Myr and 0.31% ± 0.03% at 750 
Myr, still well below the triple frequency detected ob- 
servationally in these two clusters. Thus at early times, 
the model does not form enough triples dynamically to 
match these observations. 

Furthermore the triple frequency in the Galactic field 
is observed to be significantly higher than that found 
for the dynamically formed t riples in the NGC 188 
model. Ragh avan et al.l (|2010l) find a triple frequency 
of 9% ± 2% in their complete sample of solar-type dwarf 
and subdwarf Galactic field stars. This survey combined 
spectroscopic and astrometric techniques, and is there- 
fore not limited by period; although we note that a frac- 
tion of these triples have outer periods beyond the hard- 
soft boundary in NGC 188. The maximum triple fre- 
quency reached by any of the twenty simulations in the 
N GC 188 model is 4.5% (at 7.2 Gyr). 

iTokovinin et al.l (|2006|) performed a comprehensive 
survey of short-period (1 < P < 30 days) solar- type 
Galactic field binaries within 100 pc, and find a tertiary 
{Ntripie/ Nbi nar y) frequency of 65% ± 5%. (Note that the 
denominator in this definition of the tertiary frequency is 
different from that of the triple frequency shown in Fig- 
urc !17 | and discu s sed ab ove.) Accounting for triples in the 
ITokovinin et al.l (|2006l ) sample that have outer periods 
beyond the hard-soft boundary of 10 6 5 days reduces this 
observed frequency to ^40%. However, the maximum 
tertiary frequency reached in the model for this mass and 
period range is only 10.9% ± 1.7% (at ~7 Gyr), sign ifi- 
c antly lower than observ ed by ITokovini n et all ([2006) . 

Toko vinin et al.l ()2006l ) also find the field tertiary fre- 
quency to be a strong function of the period of the in- 
ner binary, ranging from 96% (1 < P[days]< 3) to 34% 
(12 < P[days]< 30). For comparison, we find tertiary 
frequencies of 9% ± 2% and 10% ± 3%, respectively, 
at 7 Gyr for these period bins. Thus we do not find the 
same increase in the tertiary frequency for shorter-period 
inner binaries in our model. 

In conclusion, the frequency of dynamically formed 
triples produced within the NGC 188 model is signifi- 
cantly lower than triple frequencies observed in the field 
and open clusters, suggesting that real open clusters may 
form with a significant population of primordial triples. 

8.2. Period Distributions of Dynamically Formed 
Triples 

In Figure [18] we compare the orbital periods for the 
sol ar-type triples in the N GC 188 model to those of 
the ITokovinin et all (|2006l) study ( top) and to tho se in 
the Multiple Star Calatogue (MSC; lTokovininlll997[ bot- 
tom). We limit both observational samples to include 
only triples with solar-type primary stars and outer pe- 
riods of log(P2 [days]) < 6.5 (the hard-soft boundary in 
the model at the ag e of NG C 188). For comparison with 
the ITokovinin et ail (|2006l ) sample, we only include solar- 
type triples from the model with 1 < Pi(days)< 30. The 
MSC contains triple systems of all Pi, and therefore we 
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Fig. 18. — Cumulative distributions of inner (left) and outer ( right) orbital per i ods fo r solar-type dynamically formed triples in the NGC 
188 m odel compared t o the observed field triple population of Tokovinin ct al. (2006) (top) and triples in the Multiple Star Catalogue 
(MSC; Tokovinin 1997, bottom). In all plots, the line shows the distributions from the model for triples containing solar-type (0.7 Mq - 
1.1 Mq) primaries present between 6 and 7.5 Gyr . In the top panels, w e limit the model sample to include only triples with inner-binary 
periods of Pi < 30 days (to correspond with the Tokovinin et al. 2006 study). The black points in all plots show the observed samples 
including only triples with primary masses between 0.7 Mq and 1.1 Mq, and outer periods of log(P2 [days]) < 6.5 (the hard-soft boundary 
in the model between 6 and 7.5 Gyr). Both the inner and outer orbital period distributions from the model are indistinguishable from 
those of the observed samples, respectively. 



do not limit the solar-type triple sample from the model 
for this comparison. 

We note that the iTokovinin et al.l (|2006fl sample is 
nearly complete. However, the MSC is a compilation 
of triples observed in various ste llar environments from 
the literature (including those in ITokovinin et aTi r2006). 
and therefore is not meant to be a complete sample. Still 
this is the largest sample of triples in the literature, and 
therefore offers a useful comparison. 

Remarkably, both the inner and outer orbital period 
distributions for dynamically formed triples in the model 
are indistinguishable f rom those of the similar triples in 
ITokovinin et al.l (|2006l ) and the MSC, respectively. Fur- 
ther observations are needed in order define the distribu- 
tions of other orbital parameters and component masses 
for a more in-depth comparison against iV-body simula- 
tions. 

9. DISCUSSION 

The NGC 188 model matches the observed mass, 
surface-density profile, MS binary frequency and overall 
distributions of MS binary orbital elements of NGC 188 
very closely. We discuss in Section 19.11 the implications 
of this correspondence with NGC 188, given that our ini- 
tial binary population was defined empirically based on 
our observations of M35. On the other hand, the model 



creates a population of long-period circular MS- WD bi- 
naries not seen in NGC 188, and also creates too few BSs. 
Both of these highlight areas where the model requires 
improvement. We discuss these issues in Sections | 
and 19.31 respectively. 



Dynamical Evolution of the NGC 188 Solar-Type 
Binaries 



9.1. 



In Section 16.31 we found that minimal changes to the 
solar-type MS binary population occurred in our NGC 
188 model as a result of stellar dynamics, even after 7 
Gyr of evolution within the star cluster. For the solar- 
type MS binary population, dynamical encounters only 
significantly affect the long-period binaries. Most dra- 
matically, strong dynamical encounters truncate the pe- 
riod distribution at the hard-soft boundary (Figure [TT|) . 
In FigureHHwe also show that for binaries with P > 1000 
days, stellar encounters begin to significantly shift the 
eccentricity distribution to higher eccentricities. Stellar 
dynamical and two-body relaxation processes also result 
in a modest increase in the global solar-type hard-binary 
frequency, and the central concentration of binaries rel- 
ative to single stars. However for any single simulation 
(e.g, any one realization of NGC 188) no distributions of 
binary orbital parameters or masses, expect for period, 
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for the solar-type MS binaries at 7 Gyr show statistically 
significant changes as a result of dynamical processes. 

These results underscore the importance of defining a 
realistic initial binary population in star cluster simula- 
tions, as the properties of the initial hard-binary popu- 
lation persist throughout man y Gyr of dynamical evolu - 
tion (as was also discussed in iGeller fe Mathieul [2012). 
We suggest that this is best accomplished by comparison 
to observed binary populations. Here we use the solar- 
type hard binaries from the young open cluster M35 as 
a guide, although all indications to date suggest that the 
field binaries are equally appropriate. 

The NGC 188 model nearly reproduces the observed 
NGC 188 MS solar-type hard-binary population at 7 
Gyr. This result is due to the lack of dynamical pro- 
cessing of the binaries, despite having evolved within a 
fairly massive open cluster for 7 Gyr. The correspon- 
dence between the 7 Gyr MS solar-type hard binaries in 
the model and observations is because NGC 188 and M35 
themselves have very similar MS binary populations in 
both frequency and distributions of orbital parameters. 

Our results suggest that, at the age of ^180 Myr, 
NGC 188 likely had a very similar MS binary popula- 
tion as M35, and indeed both M35 and NGC 188 may 
have formed with very similar binary populations. The 
correspondence between the NGC 188 and M35 binaries 
with similar binaries in the field (which likely formed in a 
variety of different environments) introduces an intrigu- 
ing possibility that these observed binary properties may 
be ubiquitous, at least for solar-type stars. Observations 
of additional young open clusters, like M35, are essential 
to determine truly how common these binary properties 
are and the degree of scatter in the solar-type binary fre- 
quency and distributions of orbital parameters resulting 
from binary formation in different environments. Given 
the results of the NGC 188 model, we suggest here that 
even observations of older open clusters can provide valu- 
able information about the properties of young binary 
populations. 

9.2. Long-Period Circular Solar-Type-Main-Sequence - 
White-Dwarf Binaries 

In Section [6] we find a large population of solar-type- 
MS - WD binaries in the model that have periods well 
beyond the tidal circularization period, but with zero 
eccentricities (see Figures [91 and ITU]) . Such binaries are 
not observed in NGC 188 or any ot her solar-type bi- 
nary population in the literature ( e.g. IGeller fc Mathieu 
20l'a iMeibom fc Mathieul 120051: iDuauennov fc Mayor 
1991 ; iRaghavan et al.ll2010D . Essentially all of these bi- 
naries in the model started with non-zero eccentricities, 
and were later circularized during a CE episode that be- 
gan when the primary ascended the RGB or asymptotic- 
giant branch . Sim ilar binaries were also seen in the 
I Hurley et al.l ()2005l ) si mulation, and were noted as an un- 
realistic population by IGeller fc Mathieul (|2012l) in their 
comparison with NGC 188. Here we examine two general 
questions in the hopes of resolving this discrepancy with 
observations: (1) Should these binaries truly go through 
CE evolution? And, (2) if they do go through CE evolu- 
tion, should the products be MS - WD binaries in circular 
orbits with periods beyond the circularization period? 

We begin by examining whether these binaries should 
truly go through CE evolution in reality. In NB0DY6, 



as in many binary evolution codes, the decision about 
whether to send a binary with one member (or both) 
overfilling its Roche Lobe through dynamically unsta- 
ble mass transfer (e.g., CE) of thermal/nuclear mass 
transfer is controlled by the critical mass ratio q c . Gen- 
erally, in simulations utilizing q c , a quantity q\ = 
Mdonor/M accretor is calculated, and binaries with qi > 
q c undergo CE evolution, while those with qi < q c 
undergo stable mass transfer. Here we examine four 
different methods for c alculating q c from th e liter - 
ature, name l y tho s e of IHiellming fc Webbinkl (|1987|) . 
iHurlev et al.l ([20021), IChen fc Hani (120081) and the Star- 
Track code fromlBelczvnski et"aT1 (|2008l ). 

The IHurlev et al.l (120021) q„ values are the lowest 
for any given bina ry in this sample, followed by the 
IChen fc Hani ([20081) criteria for conservati v e mas s trans- 
fer, and then the IHiellming fc Webbinkl (119871) values 
used in our NGC 188 model. Thus these models 
would send all of thes e binaries through CE evolution. 
iBelczvnski et al.1 {2008) use a slightly different procedure 
to determine when CE evolution will occur, which al- 
lows for a delayed dynamical instability for such binaries 
above a critical mass ratio qadi — 3 (defined similarly to 
q c , above). If we evolve these binaries with Star Track, 
~37% avoid CE, and ^7% undergo stable mass transfer 
to become BSs. 

However, these q c values (including those used in the 
./V-body model) are highly uncert ain, and depend on how 
conservative the mass transfer is ([Chen fc Hanll2008l ). In 
the iV-body code, conservative mass transfer is assumed 
when calculating the exponents for the boundaries be- 
tween the dynamical, thermal and nuclear mass-transfer 
regimes (although non-conservative mass transfer is pos- 
sible once the regime is determined). We do not account 
for the effects of non-conservative mass transfer explic itly 
in the calculation of q c . However. I Woods et al.l (j2012f ) ar- 
gue that especially the initial stage of mass transfer from 
a giant is likely to be non -conservative due to high mass- 
loss rates. Furthermore, IChen fc Hani (|2008l ) show that 
non-conservative mass transfer effectively increases the q c 
value, which could allow additional binaries to avoid the 
CE phase, and potentially undergo stable mass transfer 
to produce BSs. 

Also, enhanced wind mass loss during the primary's 
giant phase(s) will reduce the primary mass, and may 
also increase the secondary mass through wind accre- 
tion, both of which will reduce the q\ value and poten- 
tially allow for stable mass transfer in some systems. In 
the iV-body code, we use a Reimers mass loss coefficient 
of 0.5. However the true R eimers coefficient for these 
stars may be as high a s ~1.4 (jKudritzki fc Reimerslll978l : 
ISackmann et al.lll993l) . which would allow additional bi- 
naries to avoid the CE phase, and some to become BSs. 

There is also growing evidence that this critical 
mass-ratio parameterization may be t oo simplistic 
([Woods fc Ivanovall201lUPassv et al.l l2011[). Specifically, 
IWoods fc Ivanovan^Ollh find that no parameterization 
based only on q can reliably determine the onset of unsta- 
ble mass transfer. They also find that even for systems 
that will eventually undergo CE evolution, 5% - 10% of 
the giant donor's mass may be transferred to the com- 
panion prior to this phase. For many of the systems we 
discuss here, this amount of mass would be sufficient to 
produce a BS. 
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Furthermore, Egglctonj (|2006T ) shows from an obser- 
vational sample of detached binaries possibly related to 
CE evolution that only those with progenitor mass ratios 
> 4 show the significant orbital shrinkage expected from 
a CE phase. This may suggest that only systems with 
large mass ratios undergo typical CE evolution. All the 
solar-type-MS - WD binaries from the NGC 188 model 
discussed here have qi < 2.5. Therefore perhaps most or 
even all of these binaries should avoid CE evolution in 
reality. 

In summary, these unrealistic binaries produced in the 
iV-body model may result from having incorrectly ap- 
plied CE evolution to their progenitors. In this scenario, 
these binaries could instead be BS-WD binaries at 7 Gyr 
with very similar orbital parameters and masses to the 
majority of the BSs observed in NGC 188. 

Next we investigate the second question posed above, 
namely if CE evolution is truly expected for these bi- 
naries, should the products be MS- WD binaries in long- 
period circular orbits. The post-CE orbital periods in the 
model are dete rmined by the unc ertain cxce efficiency 
parameter (seelHurlev et al.ll2002|) . In our model (as in 
IHurlev et al.l 12005) we attempt to mimic the treatment 
of llben fc Livid ( |1993f ) by setting ccce = 3. Doing so 



in the NB0DY6 treatmen t is approxim a tely e quivalent to 
using acE = 1 hi the llben fc Livid (|1993I ) treatment, 
as the two approaches use different equations for the or- 
bital and binding energies. These unrealistic MS- WD 
binaries may indicate that the acE value chosen here is 
incorrect. Perhaps the acE value (or values) for these 
binaries should result in final periods that are less than 
Pare or even mergers (potentially producing BSs) after 
the CE phase. 

Finally, in the TV-body model, we assume that all post- 
CE binaries are in circular orbits. However hydrody- 
namic simulations show that the envelope ejection can by 
asymmetric, w hich could impart a non-zero eccentricity 
to the system (Pas sv et alj|20 12). Additionally, material 
from the envelope may remain bound to t he binary after 
the ejection phase (e.g 
1201 It ISandquist et alJ 



Passv et al.l[2012tlKashi fc Sokerl 
19981 ). This material may set- 



tle into a circumbinary disk, and such post-CE bina- 
ries c an gain eccentricity through disk-binary inte rac- 
tions (jKashi fc Sokenl2011t lArtvmowicz et al.lll991[ ). In 
this scenario the product of the CE evolution may still 
be a long-period MS- WD binary, but with non-zero ec- 
centricity. 

To summarize, we identify three potential resolutions 
to this discrepancy with observations, which can be 
tested with future detailed simulations of these specific 
binaries. First, perhaps the current MS members of these 
binaries should have gained sufficient mass to become 
BSs by either undergoing stable mass transfer rather 
than CE (possibly by modifying the q c values) or by gain- 
ing mass prior to the CE phase. Second, the o>ce value in 
the model may be incorrect and perhaps should be mod- 
ified to shrink the orbits of these binaries shorter than 
Pare or even to result in mergers (possibly producing 
BSs) after the CE phase. Third, if these binaries should 
indeed undergo CE evolution and produce long-period 
MS- WD binaries, perhaps they should have non-zero ec- 
centricity. 

In closing, we note that there are roughly 9 such long- 
period circular binaries per simulation present at 7 Gyr. 



If all were indeed supposed to become BSs through mass 
transfer this would more than double the 7 Gyr BS popu- 
lation in the NGC 188 model, and importantly, these BSs 
would have similar binary characteristics to the majority 
of the NGC 188 BSs. Thus it is intriguing to consider 
the possibility that both the paucity of BSs and the over- 
abundance of long-period circular solar-type-MS - WD 
binaries may be tied to our incomplete understanding of 
the same physical mechanism. These binaries may pro- 
vide important test cases for future detailed binary evo- 
lution simulations aimed at investigating CE evolution 
and mass-transfer processes. 

9.3. The Paucity of Blue Stragglers 

As discussed in Section the NGC 188 model con- 
tains on average ~6 BSs at the age of NGC 188 with a 
detectable binary frequency of ~14.5% (P < 10 4 days). 
These values are both significantly lower than those ob- 
served in NGC 188, where we find 21 BSs with a binary 
frequency of 76% ± 19% (within the same period range). 
Here we identify potential reasons for the paucity of BSs 
in the NGC 188 model, and suggest future work that 
may further illuminate how BSs form in open clusters. 

9.3.1. Low Mass Transfer Efficiency 

iGeller fc Mathieul (|2011f ) find that the binary proper- 
ties of the long-period NGC 188 BSs (and particularly 
the companion-mass distribution) are most closely con- 
sistent with an origin in mass transfer. We find above 
that only the mass-transfer mechanism can reproduce the 
very large hard-binary frequency and distinctive distribu- 
tions of orbital periods and companion masses observed 
for the NGC 188 BSs. 

However, the mass-transfer mechanism is the least ef- 
ficient at producing BSs in the NGC 188 model. If 
we evolve the primordial binari es from our twenty sim- 
ulations in isolation (using the IHurlev et alJ 120021 BSE 
code), we would expect that at 7 Gyr on average ^5 BSs 
per simulation would have formed through mass trans- 
fer (not accounting for the expected mass loss from the 
cluster due to tidal stripping and dynamical ejections). 
Therefore the efficiency of the mass-transfer formation 
mechanism is far too low in the TV-body model than 
would be necessary to produce the ~13 BSs in NGC 188 
that we expect are the result of mass-transfer processes. 

In fact, these five mass-transfer BSs that results from 
evolving the binary population from the model in iso- 
lation are lost in a dynamical environment. Dynamical 
encounters disrupted, modified and/or ejected nearly all 
of these "proto-BS" binaries in the NGC 188 model. On 
the other hand dynamical encounters do not efficiently 
create BSs from mass-transfer processes within the NGC 
188 model. Only on average ~0.4 additional BSs per 
simulation were formed via mass transfer that were not 
expected from the population synthesis analysis. 

In Scction[7]we note that the NGC 188 model lacks the 
bimodal BS radial distribution of the true cluster, and 
instead the BSs in the model are centrally concentrated, 
with no significant halo population. This discrepancy be- 
tween the model and observations may be tied to the in- 
efficient production of BSs through mass transfer, where, 
in the model, insufficient BSs are produced in the halo to 
yield a bimodal distribution. Additional iV-body models 
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Fig. 19. — Mean number of blue stragglers between 6 and 7.5 Gyr 
as a function of the number of primordial triples in the iV-body sim- 
ulations. The values at zero primordial triples are from the NGC 
188 model. All other values are from the additional simulations 
described in Section 19.3.21 The points connected by the solid line 
show the mean total number of BSs. The points connected by 
the dashed line show the mean total number of hard-binary BSs 
(P < 10 4 days), and the points connected by the dotted line show 
the mean total number of single BSs. We find here that the total 
number of BSs increases only slightly with an increased number of 
primordial triples, but the hard-binary frequency among the BSs 
increases dramatically. Still the total number of blue stragglers in 
these simulations are a significantly lower than the number (21) 
observed in NGC 188. 

are necessary to fully understand the implications of this 
discrepancy between observations and simulations. 

We suggest that the paucity of BSs in the NGC 188 
model reflects upon an incomplete description of the 
mass-transfer mechanism in the model. The presence of 
a population of unrealistic solar-type-MS - WD binaries 
with long-periods and circular orbits further highlights 
our incomplete understanding of the interface between 
CE evolution and stable mass transfer. As discussed in 
Section 19.21 if these spurious MS- WD binaries were in- 
stead meant to undergo stable mass transfer to become 
BSs, the NGC 188 model would roughly match the ob- 
served BS number and binary properties. Our analysis 
of the BSs produ ced in the NGC 1 8 8 mo del combined 
with the results of iGeller fc Mathieu ([201 lh suggest that 
the efficiency of stable mass transfer may be significantly 
underestimated in the TV-body model (and other codes 
that use similar binary evolution prescriptions). 

9.3.2. Effects of Primordial Triples on Blue Straggler 
Production 

Another hypothesis is that the lack of BSs in the model 
is the result of a missing population of hierarchical triples 
that could produce BSs either through mergers of the in- 
ner binaries as a result of Kozai cycles a nd tidal friction 
(|Ivanoval[200l iPerets fc Fabrvckvll2009h . known as the 
KCTF mechanis m, or through collisions during dynam- 
ical encounters ([Mathieu & : Gelled l2009t iLeigh fe Sills! 
1201 lh . Kozai cycles and tidal friction are i nclud ed in the 
./V-body model as in lMardling fe Aarsethl (|2001h . Stellar 
encounters involving triples that may result in collisions 
are also modeled in detail. 

We did not include triples in the primordial population 
of the NGC 188 model, but triples do form dynamically 
(see Section [8j . In Section [7] we find that many of the 
BSs in the NGC 188 model that formed through collisions 
were members of dynamically formed hierarchical triples 
prior to the encounter that resulted in the collision. This 



is fully consistent with the analytic work of lLeigh fe Sills! 
(2011), who show that dynamical encounters involving 
triples may be common in clusters like NGC 188, and 
can lead to BS formation through collisions. 

We also investigate the NGC 188 model for dynami- 
cally formed triples that have appropriate orbital param- 
eters to potentially form BSs through the KCTF mecha- 
nism. Specifically, we search for triples with inner bina- 
ries that have a combined mass >1.15 M Q (roughly the 
lower mass li mit of BSs in NGC 188), inner periods of 
< 6 days (as in lPerets fc Fabrvcky||2009l) ■ and outer peri- 
ods between 700 and 3000 days (roughly the period range 
of the long-period NGC 188 BSs). On average, we find 
~1 such triple per simulation at any given time, and no 
BSs present in our model between 6 and 7.5 Gyr formed 
through the KCTF mechanism. Only three BSs present 
at other times within our twenty simulations formed as 
a result of Kozai cycles. 

However, it may be common for short-period (P < 

6 d ays) binaries t o fo rm with tertiary compan - 

ions (|Tokovinin et alj|2006t iFabrvckv fc Tremainell2007h , 
which are conducive to forming BSs through the KCTF 
mechanism and collisions resulting from dynamical en- 
counters. Furthermore, in Section [51 we find that the 
triple frequency in the NGC 188 model is lower than ex- 
pected from observations of real open clusters (and the 
Galactic field) , suggesting that perhaps triples should be 
included in the initial population. 

As a first step toward investigating the importance of 
primordial triples for BS production, we generated a few 
additional simulations using the same initial conditions 
as in the NGC 188 model, but also including up to 200 
primordial triples. For simplicity and to emphasize the 
types of triples that may be important for producing BSs 
in a cluster like NGC 188 through the KCTF mechanism, 
in these simulations we choose the primaries in all ini- 
tial triples to have masses of 1 M , the secondaries to 
have masses chosen randomly from a uniform distribu- 
tion between 0.2 M and 1 M , and the tertiaries to have 
masses chosen randomly from a uniform distribution be- 
tween 0.5 M^ and 2.0 M^. We choos e the inner periods 
from the iDuauennov fc Mavorl ([l991) log-normal distri- 
bution (consistent with the M35 solar- type binaries) , but 
limited to be between 2 and 50 days. The cut in inner 
period at 50 days is motivated by the obser vations of 
Galactic field triples by iTokovinin et al.l ([20061 ) who find 
that the tertiary frequency falls off quickly with increas- 
ing inner binary period. We impose a lower limit on the 
inner period of 2 days so that we don't introduce bina- 
ries that will quickly merge on their own through Case 
A mass transfer (which occurs at periods < 1.7 days for 
solar-type stars in BSE). The outer periods are chosen 
randomly between 700 and 3000 days, to reflect the pe- 
riods of the observed NGC 188 single-lined (SB1) BS bi- 
naries. The eccentricities for both orbits are chosen from 
the M35 Gaussian distribution. All angles are chosen 
randomly. 

For these simulations, we modify the initial frequency 
and period distribution of the regular (non-triple) binary 
population so that when the two components of the triple 
systems are combined with the regular binary popula- 
tion, the simulations still begin with the same log-normal 
period distribution and the initial binary frequency used 
in the NGC 188 model. There are ~200 solar-type MS 
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binaries in the initial population of the NGC 188 model 
with 2 < P(days) < 50, which sets the maximum number 
of triples that we include in these simulations. 

We ran simulations with 50, 100 and 200 primordial 
triples, and compare the results to our NGC 188 model. 
In Figure [19] we show the mean number of BSs present in 
these simulations between 6 and 7.5 Gyr. The points at 
zero triples come from our "standard" NGC 188 model 
without primordial triples. The points connected by the 
solid line show the mean total number of BSs in these 
simulations. There is a slight increase in the number of 
BSs with increasing number of triples, but even in the 
simulation containing 200 primordial triples (essentially 
replacing every solar- type binary with 2 < P[days]< 50 
in the standard NGC 188 model with a triple) the sim- 
ulation still only creates about one-third of the number 
of BSs that are observed in NGC 188. 

The BS hard-binary frequency (P < 10 4 days) in- 
creases dramatically with increasing number of triples, 
as more BSs are created within triples, and therefore re- 
tain companions. Also note that this increase in BS hard- 
binary frequency is in part due to setting the initial outer 
orbits to have periods between 700 and 3000 days, which 
are not easily disrupted. About half of the BSs in bina- 
ries that are present in these simulations between 6 and 
7.5 Gyr have origins from primordial triples. The major- 
ity of these BSs that began within primordial triples were 
involved in dynamical encounters that perturbed the sys- 
tem and in many cases exchanged in other stars prior to a 
collision that formed the BSs. Only about 25% of these 
formed in relative isolation through Kozai-induced BS 
formation. Thus, triples do play a significant role in cre- 
ating BSs through stella r co llisions (as was sug gested by 
IMathieu fc Gellerl[2009l and lLeigh fc1^l2lmh . but the 
KCTF mechanism is not a significant contributor to the 
BS population here. 

The NB0DY6 code uses the iMardling fc Aarsethl (|2001l ) 
analytic method to model Kozai oscillations and tidal 
processes in triple stars, rather than direct integra- 
tion. This method is used because direct integration 
of triple (and higher-order) stars for many Myr within 
a cluster simulation is computationally very challeng- 
ing and could also introduce further uncertainties due 
to the many necessary integration steps, each with its 
own computational u ncerta inty. It is possible that the 
IMardling fc Aarsethl (|2001f) method underestimates the 
efficiency of KCTF BS formation, which would affect the 
result s shown here. Further invest igation into the efficacy 
of the IMardling fc Aarsethl (|200l method for KCTF BS 
formation in TV-body simulations is desirable, but beyond 
the scope of this project. Additionally, further TV-body 
simulations with more realistic primordial triple popula- 
tions are required to fully test this hypothesis. 

In summary, our results here suggest that primordial 
triples may not provide the additional BSs that are re- 
quired to bring the number of BSs in the NGC 188 model 
into agreement with those observed in the true cluster. 
Specifically, even replacing all solar-type binaries with 
2 < P(days)< 50 in the standard NGC 188 model with 
triples does not reproduce the number of BSs observed 
in NGC 188. The majority of the BSs produced in these 
models through triple-mediat ed mechanisms we r e pro - 
duced by collisions. However iGeller fc Mathieu] (|2011l ) 
rule out collisions as the dominant formation mechanism 



for the majority of the NGC 188 BSs at high confidence 
(based primarily on the observed secondary-mass and or- 
bital eccentricity distributions of the NGC 188 BSs). 

9.3.3. Dependence of Blue Straggler Production on the 
Initial Binary Population and Cluster Structure 

The BS production rate is sensitive to the initial bi- 
nary parameters, and certain modifications can increase 
the number of model BSs to be consistent with the 21 ob- 
served in NGC 188. We describe the results of additional 
TV-body simulation with modified initial conditions be- 
low. However none can reproduce the observed binary 
properties of the NGC 188 BSs. 

Adding short-period binaries (e.g., by using an initial 
period distribution that is flat in log period) is the most 
efficient method for increasing the number of BSs. Here 
the BS production rate through m ergers and co l lisions 
are increased, as was also seen in the lHurlev et~aTl ([2005) 
simulation. However as discussed above such BSs do 
not match the binary properties of the NGC 188 BSs. 
Adding more high-eccentricity binaries (e.g. by using an 
initially thermal eccentricity distribution) increases the 
BS production rate through collisions, as small pertur- 
bations to a highly eccentric orbit can cause the compo- 
nents of the binary to collide. Dramatically increasing 
the overall hard-binary frequency can also increase the 
BS production rate through all mechanisms. However, 
modifying the initial period or eccentricity distributions 
or the binary frequency in this manner results in a 7 
Gyr MS binary population that is inconsistent with that 
o bserved for NGC 1 88. 

IGeller fc Malhieul (goij) find that the NGC 188 MS bi- 
naries have a modest abundance of "twin" binaries which 
is not reproduced in our model. Adding a population of 
twins to the initial binary population does not signifi- 
cantly increase the BS production rate. 

We also ran additional TV-body simulations to investi- 
gate the effect of changing the initial half-mass radius of 
the cluster (between 4 and 7 pc) and find no significant 
effect on the BS production rate at the age of NGC 188. 
After 3-5 Gyr of evolution, the clusters relax to similar 
core radii regardless of the initial values, and by 7 Gyr 
the core (and tidal) radii in all of these additional simu- 
lations are consistent with that of the NGC 188 model. 
Since the typical BS lifetime in the model is ~1.6 Gyr, 
the difference in cluster structure at a young age does 
not significantly affect the BS population at the age of 
NGC 188. 

In summary, we are unable to modify the initial con- 
ditions of the binary population or the cluster structure 
in such a way as to both reproduce the observed BS and 
MS binary populations. 

9.3.4. Blue Straggler Lifetimes 

It is possible that the paucity of BSs in the model are a 
result of underestimated BS lifetimes. Numerous uncer- 
tainties, including, for example, how efficiently fresh Hy- 
drogen can be mixed into the core, translate into uncer- 
tainties on the BS lifetimes in the NGC 188 model. For 
BSs present at 7 Gyr in the NGC 188 model, the mean BS 
lifetime is ^1650 Myr (not including prior evolution as a 
normal MS star). In order to match the number of BSs 
observed in NGC 188, we would have to approximately 
triple the lifetime of all BSs. However the total lifetime 
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of a normal MS star with a mass of a typical 7 Gyr BS in 
the model (^1.46 M©), is less than twice that of the mean 
BS lifetime. Therefore tri pling all BS lifetimes w ould be 
unphysical. Furthermore. iGlebbeek et al.l (|2008f ) suggest 
that BSs, and particularly those formed by collisions, 
will be out of thermal equilibrium and poorly mixed, and 
therefore will have shorter lifetimes than predicted by the 
equilibrium models used in the iV-body code. Finally, 
simply increasing the lifetime of all BSs would likely not 
change the BS binary frequency. 

10. SUMMARY AND CONCLUSIONS 

In this paper we present a sophisticated A^-body model 
of the old (7 Gyr) open cluster NGC 188 that matches 
the observed MS and RGB binary populations in detail, 
and thereby allows us to study the formation rate and 
mechanisms of the BS population within a realistic theo- 
retical framework. We take great care to define the initial 
conditions for the model empirically, where possible. Im- 
portantly, we employ our observations of the solar-type 
binary population of the young open cluster M35 (180 
Myr) to guide our choices for the initial binary frequency 
and distributions of period and eccentricity (and we note 
that these M35 binary characteristics are consistent with 
those of solar- type binaries in the Galactic field). 

The power of using detailed observations to both guide 
the choices for initial conditions and test the outcomes 
of A^-body open cluster simulations is evident in the ac- 
curacy with which we reproduce the observations of the 
NGC 188 solar-type binaries with our model. The MS 
binary frequency and distributions of orbital parame- 
ters agree in detail with tho se observed in the real clus- 
ter (|Geller fc Mathieul 120121) . Additionally, at the age 
of NGC 188 the model matches the observed cluster 
mass as well as the core and tidal radi i of the cluster 
dBonatto et al.l 120051: iGeller et all 120081: IChumak et "all 
12010ft . 

The MS binaries show very little evidence for modifi- 
cations by stellar dynamical encounters even after 7 Gyr 
of evolution within the cluster environment. Indeed only 
the long-period binaries show evidence for orbit mod- 
ifications from stellar encounters. Specifically, strong 
dynamical interactions are responsible for breaking up 
the initially very long-period binaries, and establishing 
a hard-soft boundary at P ~ 10 6 5 days. Stellar en- 
counter (most often fly-by encounters) begin to shift the 
eccentricity distribution for MS-MS binaries to higher 
eccentricities for binaries with P > 1000 days, and this 
effect increases towards longer periods. Two-body re- 
laxation processes also result in a ^20% increase in the 
detectable binary frequency over the 7 Gyr of evolution 
and the central concentration of the MS binaries in the 
cluster at the age of NGC 188. However, in general, the 
short-period MS binaries in the NGC 188 model show 
little signs of dynamical processing, which further em- 
phasizes the importance of defining an accurate initial 
binary population for constructing realistic open cluster 
models, as most of the hard-binary characteristics persist 
throughout the evolution of the cluster. 

This finding may also suggest that observations of the 
present-day binaries in open clusters can provide valuable 
information on their primordial binary populations. The 
similarities between the solar-type binary populations in 
NGC 188 (at 7 Gyr), M35 (at 180 Myr) and the Galac- 



tic field are remarkable, and may indicate that binaries 
form with similar characteristics within a variety of en- 
vironments. Comparable studies of binary populations 
in additional rich open clusters are necessary to explore 
this further. Nevertheless, we suggest that the observa- 
tionally defined distributions of binary orbital parame- 
ters and binary frequency used here should be strongly 
considered for use in future open cluster models. 

However, despite having matched the observed MS bi- 
naries, the NGC 188 model produces only one-third as 
many BSs, and these BSs have only one-fifth the hard- 
binary frequency (P < 10 4 days), as compared to the BSs 
observed in NGC 188. For the full 7 Gyr of the modeled 
cluster, mergers and collisions dominate BS formation, 
with collisions producing the majority of the BSs at the 
age of NG C 188. This is in stark c ontrast to the obser- 
vations of Ge ller fc Mathieul ()2011f ) who show that the 
secondary-mass distribution of the majority of the BSs 
in NGC 188 is inconsistent with the collisionally pro- 
duced BSs in the model. Rather the binary properties of 
the NGC 188 BSs, and particularly the secondary-mass 
distribution, point to a mass-transfer origin. 

A potential explanation for the missing BSs may also 
resolve a second discrepancy between the model and 
observations, namely the population of circular MS- 
WD binaries with periods well beyond the tidal circu- 
larization period created within the model as a result 
of CE evolution. Such bi naries are not observed in 
NGC 188, the G alactic field dDuauennov fc Mavonll99l"l : 
Ragh avan et al.l 12010ft or the ma ny oth er open cluster 
surveyed by iMeibom fc Mathieul (|2005l ). Perhaps the 
MS stars in these binaries were instead meant to become 
BSs by entering a phase of stable mass transfer. This is 
a compelling solution as the number of unphysical MS- 
WD binaries in the model is consistent with the number 
of missing BSs, and the binary characteristics of these 
BSs would be consistent with the majority of the NGC 
188 BSs. 

The creation of these unphysical long-period circular 
MS- WD binaries in the model highlights our incomplete 
understanding of CE evolution, and particularly which 
binaries should enter a CE phase, which should enter 
stable mass transfer and what the products of CE evo- 
lution should look like. Improvements may be necessary 
in the treatment of mass transfer and CE evolution in 
models like NB0DY6 and BSE, which may also have impor- 
tant implications for formation rates of other interacting 
binary populations, such as W UMa's, CVs, symbiotic 
stars, etc. These specific unphysical binaries created in 
the NGC 188 model may provide important test cases 
for future detailed binary evolution simulations aimed at 
studying CE evolution and mass-transfer processes. 

Missing BSs in the NGC 188 model may also 
point to a mi ssing population of pr imordial triples. 
The Pleiades M ermilliod et al.l Il992j ) and Praesepe 
(|Mermilliod fc Mavorl 19991 ) both have significantly 
higher observed triple frequencies than the NGC 188 
model at the respective cluster ages, as does the Galactic 
field (iDuquennoy fc Mavorl 119911: iTokovinin et~aT1 120061: 
iRaghavan et al.l | 2010) . However, our additional simu- 
lations (Section 19.3. 2|) suggest that adding primordial 
triples will not dramatically increase the number of BSs 
present at 7 Gyr. Still, additional A^-body simulations, 
with more realistic initial triple populations, are required 
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to fully investigate if primordial triples can contribute 
significantly to BS production. 

In closing, we note that only recently have both the ob- 
servations and iV-body simulations matured sufficiently 
to permit such detailed comparisons as we present here 
between true binaries in a dynamically evolved open clus- 
ter and analogous binaries simulated in a realistic iV- 
body open cluster model. Open clusters are unique in 
their accessibility to both detailed comprehensive obser- 
vational surveys and full-scale sophisticated -ZV-body sim- 
ulations, opening a window into the dynamical evolution 
of star clusters that is only now becoming reachable. Fur- 
thermore, advances in both hardware and software now 
readily allow for multiple large open cluster simulations 
to run simultaneously (as we have done here), which en- 
able more statistically robust analyses and a broader ex- 



ploration of parameter space ()Mathieul [2008). This ca- 
pability, combined with the upcoming results from our 
WOCS observations of a number of open clusters with 
a wide range in age, will continue to advance our un- 
derstanding of star cluster dynamics, the origins of BSs 
and the roles that binary and multiple stars play in star 
cluster evolution. 
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